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ABSTRACT 

The evolution of the 2006 outburst of the recurrent nova RS Ophiuchi was followed with 12 X- 
ray grating observations with Chandra and XMM-Newton. We present detailed spectral analyses 
using two independent approaches. From the best dataset, taken on day 13.8 after outburst, we 
reconstruct the temperature distribution and derive elemental abundances. We find evidence for at 
least two distinct temperature components on day 13.8 and a reduction of temperature with time. 
The X-ray flux decreases as a power-law, and the power-law index changes from —5/3 to —8/3 around 
day 70 after outburst. This can be explained by different decay mechanisms for the hot and cool 
components. The decay of the hot component and the decrease in temperature are consistent with 
radiative cooling, while the decay of the cool component can be explained by the expansion of the 
ejecta. We find overabundances of N and of a elements, which could either represent the composition 
of the secondary that provides the accreted material or that of the ejecta. The N overabundance 
indicates CNO-cycled material. From comparisons to abundances for the secondary taken from the 
literature, we conclude that 20-40% of the observed nitrogen could originate from the outburst. The 
overabundance of the a elements is not typical for stars of the spectral type of the secondary in the 
RS Oph system, and white dwarf material might have been mixed into the ejecta. However, no direct 
measurements of the a elements in the secondary are available, and the continuous accretion may 
have changed the observable surface composition. 

Subject headings: novae, cataclysmic variables - stars: individual (RSOph) - X-rays: stars - shock 
waves - methods: data analysis - binaries: symbiotic 
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1. INTRODUCTION 

Nova explosions occur in binary systems containing a 
white dwarf (WD) that accretes hydrogen-rich material 
from its companion. When 10~ 6 — 10 _4 M Q have been 
accreted (depending on the WD mass), ignition condi- 
tions for explosive nuclear burning are reached and a 
therm onuclear runaway (TNR) occurs (fStarrficl d et al.l 
2008). Material dredged up from below the WD surface 
is mixed with the accreted material and violently ejected. 
While nuclear burning continues, the WD is surrounded 
by a pseudo atmosphere, and the peak of the spectral 
energy distribution (SED) shifts from the optical to soft 
X-rays as the radius of the p seudo photosphere shrinks 
(jGallagher fc Starrfieldl[T978f ). Observations of novae in 
soft X-rays therefore generally yield no detections until 
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the photosphere recedes to the regions within the outflow 
that are hot enough to produce X-rays. For some novae 
this has been observed, and the X-ray spectra during 
this phase rese mble the class of Super Soft X-ray Binary 
Sources fSSS. iKahabka fc van den Heu^fl997n . This 
phase is therefore called the SSS phase. 

Observational evidence (from optical observations) and 
theoretical calculations indicate two abundance classes 
of novae, those with overabundance of carbon and oxy- 
gen (CO novae) and those with ove rabundance of oxy- 
gen and neon (ONe n ovae; see, e.g.. lAndrea et ani!994t 
I Jose fc Hernanj 119981) . Since the pressure on the white 
dwarf surface is not high enough for the production of 
C, O, or Ne during the nova outburst, these abundance 
classes reflect the composition of the WD. This indicates 
that core material is dredged-up into the accreted ma- 
terial and the gases are mixed before being eje cted into 
space (jStarrfield et allfioM iGehrz et aLlflOOSh . In ad- 
dition to dredged-up WD material, the ashes of CNO 
burnin g during the outburst have frequently be en ob- 
served (|Andrea et al.lll994t Uose fc Hernand ll998). The 
composition of the ejected material is thus highly non- 
solar. 

RS Oph is a Recurrent Symbiotic Nova, which erupts 
about every 20 years. The l atest outburst occurr ed on 
2006 February 12.83 f=davO; lHTrosawa et al.ll2006f ). The 
mass donor is a red giant (M2III), and the expanding 
ejecta interact with the pre-existing stellar wind setting 
up shock s ystems. The compositio n of the red giant was 
studied by iPavlenko et alj (|2008l ) , who found that the 
overall metallicity does not seem to be significantly dif- 
ferent from solar ([Fe/H] = 0.0 ±0.5), C is underabun- 
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dant ([C] = -0.4), and N overabundant ([N] = +0.9). 
UV spectra taken with IUE during the 19 85 outburst 
provi ded evidence that N was overabundant (|Shore et al.1 
199G). Lines of C were observed, b ut no detaile d abun - 
dance analyses were carried out by IShore et al.1 (|1996f h 
IContini et al.1 (|1995l ) determined an N/C abundance ra- 
tio of 100 and N/H=10 from optical spectra taken on 
day 201. From their absolute abundance of N and Fe, an 
abund ance rati o of N/ Fe=15 relative to solar can be de- 
rived. ISniidersI (|1987af) found N/0=1.1 and C/N=0.16, 
and they caution that the evolved seco ndary can already 
be C/N depleted. IContini et al.l (|1995h found significant 
underabundance of O/H and of Nc/H of ~ 10% solar but 
high abundance ratios of Mg/Fe=5.4 and Si/Fe=7.2. 

During the first month after outburst, intense hard 
X-ray emission, that originated from the shock, was 
observed with Swift and the Ro s si X-ray Timing Ex- 
plorer CRXTE: iBode et alj l2006t lOsborne et all 120081 : 
ISokoloski et al.ll2006h . Swift X-Ray Telescope (XRT) ob- 
ser yations carri e d out between days 3-26 were analyzed 
by IBode et all (|2006[ ) who applied single-temperature 
MEKAL models to the X-ray spectra. They deter- 
mined temperatures and wind column densities, Nyy = 
./Vh (total) — Ah (interstellar) . The interstellar value of 
./Vh (interstellar) = 2.4 x 10 21 cm~ 2 has been deter- 
mine d from Hi 21cm measurements (Hicllmin g et al.1 
1986). This value is consistent with the visual extinction 
{E{B - V) = 0.73 ± 1) dete r mined from IUE ob serva- 
tions in 1985 (|Snijderslfl987bT l. IBode et al.l (|2006h con- 
verted the temperatures found from the MEKAL models 
into derived shock velocities v s , assuming that the X- 
rays were produced in the blast wave driven into the cir- 
cumstellar material following the outburst. Before day 
~ 6 after outburst they found a power-law decay t~ a 
with an approximate index a = 0.6, 0.5, and 1.5 for w s , 
iV\v, and the flux (unabsorbed, i.e., corrected for inter- 
stellar absorption), respectively. These results compare 
well with model predictio ns of the RS O ph system pre- 
sente d by O'Brien et al. (|1992l - see also IBode fe Kahnl 
1985). According to these models, the evolution can be 
divided into three phases. The first phase (I), where 
the ejecta are still important in supplying energy to the 
shocked stellar wind of the red giant, lasts only a few 
days. The second phase (II) commences when the blast 
wave is being driven into the stellar wind and is effec- 
tively adiabatic. This phase is expected to last until the 
shocked material is well cooled by radiation (phase III). 
The physics behind these phases of evolution, together 
with the density distribution in the wind, determine the 
evolution of temperature with the corresponding veloc- 
ity of the shock, un absorbed fluxes, an d the absorbing 
c olumn of the win d (IVaytet et all I2007D . 

ISokoloski et al.l (|2006h analyzed X-ray data taken 
between days 3-21 with RXTE, and from thermal 
bremsstrahlung models found that the temperature de- 
creased with time t as t~ 2 / 3 . They concluded that the 
speed of the blast wave produced in the nova explosion 
decreased with t -1 / 3 . However, the RXTE data with 
their low sensitivity at low energies did not favor the 
measurement of the wind column density Ny/ . 

A Chandra High Energy Transmission Grating Spec- 
trograph snapshot of the blast w ave obtained at the 
end of day 13 and analyzed by iDrake et al.l (2008) 



shows asymmetric emission lines sculpted by differen- 
tial absorptio n in the circumste llar medium and explo- 
sion ejecta. IDrake et al.l ([2008) found the lines to be 
more sharply peaked than expected for a spherically- 
symmetric explosion and concluded that the blast wave 
was collimated in the direction perpendicular to the line 
of sight, as al so suggested by contemporaneous radio in- 
terferometry (jO'Brien et al.ll2006l) . 

The SSS phase was observed aft er day ~ 30 and ende d 
before day ~ 100 after outburst (|Osborne et al.ll2008h . 
Three high-resolution X-ray spectra were taken durin g 
this phase which are described by N ess et all ([2007) . 
The SSS emission longward of ~ 12 A (E > 1 keV) out- 
shines any emission produced by the shock at these wave- 
lengths, however, all emission s hortward of 12 A orig i- 
nates exclusi vely from the shock (jNess et al.ll2008L 120071 1 . 
We note that IBode et all (|2008l ) show tentative evidence 
for emission between 6-12 A that may reflect the evolu- 
tion o f the SSS. The SSS spectra analyzed by lNess et al.1 
(|2007[ ) contain emission lines on top of the bright SSS 
continuum which, combined with blue-shifted absorption 
lines, were first attributed to P Cygni profiles ( Ness et all 
2006), but may also originate from the shock iNess et alJ 
(f2007L 

An analysis of all X-r ay grating spectra was presented 
bv iNelson et al.l (|2008l ). They discovered a soft X-ray 
flare in week 4 of the evolution in which a new sys- 
tem of low-energy emission lines appeared. With their 
identifications of the emission lines, they derived veloc- 
ities of 8,000 — 10, 000 km s -1 which is consistent with 
the escape velocity of the WD, and the new component 
may thus represent the outflow. From preliminary at- 
mosphere models they also determined the abundance 
ratio of Carbon to Nitrogen of 0.001 solar. This is a 
factor 10 lower tha n C/N abundance measure ments by 
IContini et al.1 (|1995f ) and a factor 100 lower than lSnijdersI 
From He-like line flux ratios they confirm that 



the shock plasma is collisionally dominated. They mea- 
sured line shifts and line widths and found that the mag- 
nitude of the velocity shift increases for lower ionization 
states and longer wavelengths. In addition, as the wave- 
length increased, so did the broadening of the lines. They 
discuss bow shocks as a possible origin for the line emis- 
sion seen in RS Oph. From multi-t emperature plasma 
modelling of the early X-ray spectra. [Nelson et al.l ( 2008) 
needed four temperature components. While they found 
reasonably good reproduction of the Chandra spectrum, 
the same model was in poor agreement with the simulta- 
neous XMM-Newton spectrum. Their model underpre- 
dicts lines of O and N, and they concluded that these 
elements are overabundant, and that the lines originated 
in the ejecta. 

The structure of this paper is as follows: In fj2] we 
present 12 X-ray grating observations taken between 
days 13.8 and 239.2 after outburst, focusing only on the 
emission produced by the shock. We measure emission 
line fluxes and line ratios in Ej3j and in S]4]we present sup- 
porting models. We compute multi-temperature spectral 
models with the fitting program xspec f £|4.2p and recon- 
struct a continuous temperature distribution based on a 
few selected emission lines ( §4.3|) . yielding the elemental 
abundances. In ij4.4l we compare the results of these two 
model approaches. We dedicate a separate section f ^4.5|) 
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to the discussion of systematic uncertainties, as all given 
error estimates are only statistical uncertainties. In $5] 
we discuss our results and summarize our conclusions in 

m 

2. OBSERVATIONS 

In this paper we analyze five X-ray grating spec- 
tra taken with Chandra and five with XMM-Newton. 
We use the High- and Low Energy Transmission Grat- 
ing spectrometers (HETG and LETG, respectively) 
aboard Chandra and the Reflection Grating Spectrom- 
eters (RGS1 and RGS2) aboard XMM-Newton to obtain 
data between 1 A and 40 A. In Table [1] we list the start- 
and stop times, the corresponding days after outburst, 
the mission and instrumental setup, ObsIDs, and net ex- 
posure times for each observation. We hav e extracted 
the sp ectra in the same way as described by iNess et ail 
(2007) using the standard tools provided by the mission- 
specific software packages SAS (Science Analsis Software, 
version 7.0) and CIAO (Chandra Interactive Analysis of 
Observations, version 3.3.0.1). While pile up in the zero- 
th order of the Chandra HETGS observation may lead 
to problems in cen troiding the extract ion regions for the 
dispersed spectra (iNelson et al.l 120081 ). we are confident 
that the standard centroiding is accurate enough. For 
example, the wavelengths of strong lines from the two 
opposite dispersion orders agree well with each other. 

We have also extracted spectra from the XMM-Newton 
European Photon Imaging Camera (EPIC), concentrat- 
ing on the observations recorded with the Metal Oxide 
Semi-conductor chips (MOS1). We have used standard 
SAS routines for the extraction of spectra and have cor- 
rected for pile up using annular extraction regions that 
avoid extracting photons from the innermost regions of 
the point spread function (PSF) following the instruc- 
tions provided by the SAS software. We only need the 
MOS1 data for ObsID 0410180201 (26.1 days after out- 
burst) in ij4.21 for which an inner radius of 300 pix (15") 
has to be excluded to avoid significant pile up. 

2.1. Description of spectra 

An overview of all grating spectra is presented in 
Fig. [1] with the instrument and day after outburst in- 
dicated in the legends of each panel. All spectra taken 
on days 13.8 and 26.1 after outburst (top three panels) 
are characterized by a hard, broad con tinuum spectrum 
with additional s trong emission lines (|Ness "eTaTI [20061 
iDrake et al.ll2008l ). The count rate on day 26.1 is signif- 
icantly lower than that on day 13.8, and the shape of 
the continuum is different. On day 26.1 a new com- 
pone nt is observed longward of ~ 20 A ([Nelson et al.l 
I2008D that could be associated with the SSS spectrum 
that was clearly det ected three days later with Swift 
(|Osborne et alj [20081 ). However, the spectral shape of 
this new component on day 26.1 is quite different from 
the spectra observed on days 39.7, 54.0, and 66.9 (next 
three panels). These spectra are dominated by the SSS 
spectrum (jNess et al.ll2007D between 14 A and 37 A, while 
the em ission from the shock dominates shortward of 
~ 15 A (INess et al.ll2008D . After day ~ 100, the SSS spec- 
trum has disappeared, and those spectra display emission 
lines with a weak continuum. The short- wavelength lines 
are only seen in the early spectra while those between 



TABLE 1 
Grating observations of RS Oph 



Date 


Day a 


Mission 


Grating ObsID 


exp. time 


start-stop 






/ detector 


(net; ks) 


Febr. 26, 15:20 


13.81 


Chandra 


HETG 


7280 


9.9 


-Febr. 26, 18:46 


13.88 




/ACIS 






Febr. 26, 17:09 


13.88 


XMM 


RGS1 


0410180101 


23.8 


-Febr. 26, 23:48 


14.2 




RGS2 




23.8 


March 10, 23:04 


26.1 


XMM 


RGS1 


0410180201 


11.7 


-March 11, 02:21 


26.3 




RGS2 




11.7 


March 24, 12:25 


39.7 


Chandra 


LETG 


7296 


10.0 


-March 24, 15:38 


39.8 




/HRC 






April 07, 21:05 


54.0 


XMM 


RGS1 


0410180301 


9.8 


-April 08, 02:20 


54.3 




RGS2 




18.6 


April 20, 17:24 


66.9 


Chandra 


LETG 


7297 


6.5 


-April 20, 20:28 


67.0 




/HRC 






June 04, 12:06 


111.7 


Chandra 


LETG 


7298 


19.9 


-June 04, 18:08 


111.9 




/HRC 






Sept. 06, 01:59 


205.3 


XMM 


RGS1 


0410180401 


30.2 


-Sept. 06, 17:30 


205.9 




RGS2 




30.2 


Sept. 04, 10:43 


203.6 


Chandra 


LETG 


7390 


39.6 


-Sept. 04, 22:26 


204.1 




/HRC 






Sept. 07, 02:37 


206.3 


Chandra 


LETG 


7389 


39.8 


-Sept. 07, 14:29 


206.8 




/HRC 






Sept. 08, 17:58 


207.9 


Chandra 


LETG 


7403 


17.9 


-Sept. 08, 23:36 


208.2 




/HRC 






Oct. 09, 23:38 


239.2 


XMM 


RGS1 


0410180501 


48.7 


-Oct. 10, 13:18 


239.7 




RGS2 




48.7 



"after outburst (2006, Feb. 12.83) 



12 A and 25 A can be seen in all spectra, however, with 
different relative strengths. 

In Fig. [2] we show the X-ray spectrum taken on day 
13.8. For this plot we have converted the number of 
counts in each spectral bin to photon fluxes, simply di- 
viding the number of counts by the effective areas ex- 
tracted for each spectral bin from the instrument cali- 
bration. With grating spectra such a conversion is suf- 
ficiently accurate because of the precise placement of 
the recorded photons into the spectral grid. In con- 
trast to low-resolution X-ray spectra taken with CCDs, 
the photon redistribution matrix of grating spectra is 
nearly diagonal. Below 16.5 A we show the Chan- 
dra/MEG spectrum, and above this wavelength, where 
the MEG has extremely low sensitivity, the combined 
XMM-Newton/RGS spectra are shown. The strongest 
lines seen in the spectrum originate from H-like and He- 
like ions of Sxvi and Sxv (4.73 and 5.04 A), Sixiv and 
Si XIII (6.18 and 6.65A), Mgxn and Mgxi (8.42 and 
9.2 A), Nex and Neix (12.1 and 13.5 A), Ovm and O vn 
(18.97 and 21.6 A), andNvn (24.78 A). Also some of the 
3p-ls lines are detected, e.g., Mgxn at 7.11 A, Mgxi at 
7.85 A, and Ovm at 16 A. The H-like and He-like lines 
of elements with higher nuclear charge arise at shorter 
wavelengths, and strong lines at short wavelength in- 
dicate high temperatures. Several Fe lines are present, 
e.g., Fexxv (1.85+1.86+1.87 A) and Fexxiv at 10.62A 
as well as low-ionization lines of Fexvn at 15.01 A and 
12.26 A. These lines cannot be formed in the same region 
of the plasma and it is thus not isothermal. 

In Fig. [3] we show the combined XMM-Newton/RGS 
spectra taken on day 26.1, in the same units as in Fig. [3 
for direct comparison. While on day 13.8 the strongest 
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Fig. 1. — 11 X-ray grating spectra extracted from 12 Chandra and XMM-Newton observations taken on the dates and with the 
instruments given in the legends. Plotted are the raw spectra in counts per bin per individual exposure time (see Table [TJ- Bin sizes are 
0.005 A, 0.0025 A, and 0.01 A for MEG, HEG, and RGS1, RGS2, and LETG, respectively. Four Chandra spectra taken between days 204 
and 208 after outburst are combined. RGS1 and RGS2 spectra are combined. 
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Fig. 2. — X-ray spectrum of RS Oph on day 13.8, in photon flux units, taken with Chandra /MEG shortward of 16.5 A and with XMM- 
Newton/KGS longwards. We label H-like and He-like lines in italic font with the H-like lines in bold-face. Other lines are labeled with 
roman font. The high ionization stages indicate temperatures up to 10 8 K while the additional presence of low-ionization stages show that 
the plasma is not isothermal. 
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Fig. 3. — XMM-Newton/TiGSl and RGS2 spectra (combined) taken on day 26.1 and converted to photon flux units. The strongest 
emission lines are labelled as explained in Fig. [5] 

lines are formed at wavelengths shortward of 10 A, the 
Nex line at 12.1 A is now the strongest line. This could 
mean that the temperature and/or the neutral hydro- 
gen column density have decreased. The relative line 
strengths of H-like to He-like lines are significantly lower 
for all elements (see, e.g., Mgxn to Mgxi). This is 
clearly a temperature effect, and the plasma is cooling. 
Longwards of 25 A a new component can be seen. The 
fact that only three days later the SSS spectrum was 
observed with Swift (jOsborne et al.|[2008h suggests that 
this emission represents the onset o f the SSS phase (e.g., 
iBode et al.|[2006tlNelson et al.ll2Q08D . However, while the 
SSS spectra observed on day 39.7 range from ~ 15 — 30 A 



(Fig.[T]), the RGS spectra shown in Fig. [3] show only ex- 
cess emission longward of ~ 20 A (see also Fig. At 
23.5 A a deep absorption edge fr om Oi has been fo und 
in the SSS spectra of RS Oph by iNess et all (|2007f) (see 
also bottom panel of Fig. [5]). The hard portion of an 
early faint SSS spectrum might be entirely absorbed by 
circumstellar neutral oxygen in the line of sight, while 
the shock-induced emission may originate from further 
outside, thus traversing through less absorbing material. 
Also, in the standard picture of nova evolution, the peak 
of the SED is expected to shift from long wavelengths 
to short wavelengths while the radius of the photosphere 
recedes to successively hotter layers, and the observed 
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emission would be consistent with this picture. How- 
ever, the spectrum has more chara cteristics of an emis - 
sion line spectrum (see Fig. 151 and INelson et all [2008) . 
but only the lines at 24.79A and 28.78+29.1+29.54 A 
can be identified as Nvn and as the Nvi He- like triplet 
lines, respectively. In between the Nvn and Nvi lines 
no strong lines are listed in any of the atomic databases. 
The strongest emission line in this range is observed as 
a narrow line at 27.7 A (FWHM 0.08 A) with a line flux 
of (2.7 ± 0.4) x 10 _13 ergcm _2 s _1 . The only possible 
identifications would be Arxiv (27.64 A and 27.46 A) or 
Caxiv (27.77 A). Both appear rather unlikely identifi- 
cations, as no Ar lines are detected in any of the other 
spectra, and for Caxiv, stronger lines are expected at 
24.03 A, 24.09 A, and 24.13 A, but are not detected. A re- 
markable aspect is that the 27.7-A line is so narrow while 
the N vn line shows an extremely broad profile (see ^3. 1[) . 
Another unidentified line is measured at 23.6 A, but we 
experience the same difficulties in finding an identifica- 
tion. This could be residual continuum emission if the 
abso rption feature at 23. 5 A is interpreted as interstellar 
Oi. INelson et al.l (|2008h suggested that some of these 
lines are blue-shifted Nvi and Cvi lines, but this re- 
quires extremely high velocities and is in contradiction 
to the non-detection of the Cvi Lya line and the low 
C abundance reported in the same paper. In any case, 
this component is likely not part of the shock systems, 
and the discussion of it is beyond the scope of this pa- 
per. While the Ovni and Ovu lines might be part of 
the shock, we treat the interpretation of these lines with 
care. We will also include the Nvn line in our analyses 
as if it were formed in the shock, and any inconsisten- 
cies based on this line can be understood as supporting 
evidence that this component is unrelated to the shock 
emissio n. For more detail s of this component we also 
refer to INelson et all (|2008f) . 

In Fig. 2] we show the photon flux spectrum taken with 
Chandra LETGS on day 39.7. The SSS spectrum dom- 
inates all emission longward of ~ 14.5 A, and we only 
show the wavelength range relevant for this paper. The 
ratio of H-like to He-like lines is lower than in the ear- 
lier spectra, indicating that the temperature has con- 
tinued to decrease. Since the lines are formed short- 
wards of the high-energy (Wien) tail of the SSS spectrum 
(14.5 Aw 0.86 keV), they are not affected by photoexci- 
tations and originate exclusively from the shock. 

In Fig. [5]we show one of the spectra taken after the SSS 
had turned off. All emission lines are significantly weaker 
and the ratio of H-like and He-like lines is again lower 
than in the previous observation. All short-wavelength 
lines are extremely weak or are not detected. 

Next we integrate the photon flux spectra over the 
range 7-11 A (l.l-1.8keV) in order to obtain X-ray 
fluxes. We do not correct for absorption, thus yield- 
ing fluxes at Earth. Since the fluxes are extracted from 
above 1 keV, the effects from absorption are small, and 
particularly the relative evolution of the absorbed and 
non-absorbed fluxes is the same. The wavelength range 
over which the fluxes are integrated is a compromise be- 
tween collecting as much information as possible from 
the observations before day 39.7 and after day 66.9 while 
excluding as much as possible of the emission from the 
SSS on days 39.7, 54.0, and 66.9. The results are illus- 



trated as a function of time in Fig. [6l For comparison 
we include rescaled Swift/XRT count rates (0.25-10 keV) 
taken after day 106. At this late stage of the evolution, 
the spectral shape hardly changes (see Table [5]) , yield- 
ing a direct correlation between X-ray flux and count 
rate. With the assumption of no spectral changes be- 
tween days 106 and 250, we can also use the count rates 
integrated over the full Swift XRT band pass, as addi- 
tional emission in the larger wavelength range also scales 
directly with the count rate. Since we are not interested 
in the absolute flux from the Swift observations, we chose 
a scaling factor of 3 x 10 -12 erg cm -2 s _1 cps -1 to yield 
the same values as the grating fluxes for days 111.7-239.2. 
The rescaled XRT count rate follows the same trend as 
the fluxes obtained from the grating spectra. We include 
two power-law curves, and the early evolution evolves 
more like i~ 5 / 3 , while the later evolution (after day 100) 
clearly follows a t~ 8 / 3 trend. We observe the same be- 
havior if we use a larger wavelength range and exclude 
the observations between days 39.7 and 66.9. 

3. MEASUREMENT OF EMISSION LINES 
3.1. Line shifts and profiles 

In order to determine velocities from the emission lines 
we have measured wavelengths and line widths in excess 
of the instrumental line broadening function for a num- 
ber of strong lines with well-identified rest wavelengths, 
Xq. For the narrow wavelength range around the lines we 
have accounted for continuum emission by defining a con- 
stant local offset on top of the instrumental background 
that can be treated as an 'uninteresting' free parame- 
ter. For each line j we have used a normalized Gaus- 
sian profile with wavelength Xj and line width Oj and 
folded this profile through the instrumental response us- 
ing the ID L tool scrmf provi ded by the PINTofAle pack- 
age (|Kashvap fc Dr ake 2000) before comparing with the 
measured count spectra. We have determined the statis- 
tical measurement uncertainties for Xj and o~j from the 
2x2 Hesse matrix as defined in Eq. [5]of the appen dix sec- 
tion, w hich is based on an approach proposed by IStrond 
(1985). Systematic uncertainties are difficult to assess 
and are not included in our error estimates (for details 
see N4.5[) . Those can arise from fluctuations in the under- 
lying continuum and line blends. While the former has 
a stronger effect on weak lines, the latter can affect any 
line. For this reason we chose lines for which no strong 
nearby lines are known to arise. We have iterated Xj 
and o~ j , and in each iteration step we have adjusted the 
normalizatio n utilizing the fixed point iteration scheme 
described by Ness fc Wichmanr] (|2002] ). The normaliza- 
tion factor can be converted to line fluxes (see £|3.2[) . 

The results are listed in Table [2] The line shifts 
Xj — Xq and Gaussian line widths <jj (both measured 
in mA= 10 _3 A) are converted to corresponding Doppler 
velocities using the rest wavelengths Ao listed in the first 
column. The measurement uncertainties of line shifts 
and widths are correlated uncertainties, and account for 
the uncertainties in the respective other values. The Ne x 
line at 10.23 A is relatively weak in all observations, and 
the results from this line may be less certain due to ad- 
ditional systematic uncertainties from fluctuations in the 
underlying continuum. The Fexvn line could be blended 
with the weak Ovni ls-4p (A = 15.18 A) line, and the 
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Fig. 4. — Same as Fig.[3]for the Chandra/LETG spectrum taken on day 39.7. Longward of ~ 14.5 A the SSS spectrum dominates. 
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Fig. 5. — Same as Fig. [3]for the Chandra/LETG spectrum taken on day 111.7. The SSS spectrum has disappeared and emission lines 
longward of ~ 14.5 A can be observed. 



accuracy of the results from this line might suffer from 
line blending. All lines measured from observations taken 
after day 26.1 are weaker, and the uncertainties on the 
results from these observations have to be increased by 
at least 20% due to fluctuations in the continuum. 

In Fig. [7] we illustrate the measured line shifts (top 
left) and widths (top right) and the corresponding ve- 
locities (respective bottom panels) for the observations 
taken on day 13.8. All lines are significantly blue-shifted, 
the short-wavelength lines by 200 — 800 km s -1 and the 
lines of O vm, O vn, and Nvii at longe r wave lengths 
by more than 1200 km s -1 . iNelson et all (|2008f) found 
similar values and concluded that there was a trend of 
increased velocities with wavelength and thus with for- 
mation temperature. However, with a different set of 



lines we come to a different conclusion. First, we have 
not used the unresolved He triplet lines of Mgxi and 
Si XIII to avoid additional systematic uncertainties from 
line blends (see £|4.5|) . Then, we have included the Nex 
and Fc xvn lines that lie in between the Mg lines and the 
O VIII line. Although we caution that these lines may suf- 
fer from additional systematic uncertainties, there seems 
to be more of an abrupt change rather than a systematic 
trend with these additional lines included. Interestingly, 
the lines with large r blue shifts o rigina te only from oxy- 
gen and nitrogen. iDrake et al.l (|2008| ) investigated the 
possibility that the line profiles are dominated by com- 
plex absorption patterns in their red wings, leading to 
apparent blue shifts. In that case the column density in 
the respective line is a stronger driver for line shifts than 
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the temperature, and oxygen and nitrogen might exhibit 
deeper column densities than other elements, possibly 
owing to higher elemental abundances. 

The line widths are all about 800 - 1000 km s -1 , with 
the exceptions of Nex (10.23 A) and Fexvn (15.01 A 
and 16.78 A) which are narrower (bottom right panel of 
Fig. [7|). Since fluctuations in the continuum and line 
blending cannot lead to narrower lines, it is not clear 
to us why these particular lines are narrower, but we 
cannot confirm a trend with long-wavelengt h lines being 
broader as reported by iNelson et al.l (|2008| ). Nelson et 
al. seem not to have accounted for the instrumental line 
broadening when computing line widths, but the instru- 
mental line broadening is only ~ 0.01 A for the MEG 
and ~ 0.03 A in the RGS. The instrumental line profile 
is roughly Gaussian, and since the convolution of two 
Gaussians is again a Gaussian, the resulting line width is 
dominated by the broader line, and in the case of most 
lines, the instrumental line broadening can be neglected. 

After day 13.8, all line shifts except those for the 
Ovni and Nvn lines fluctuate around the same value 
of - 500 - 800 kms" 1 (see Table HJ). The O and N lines 



that show extreme blue-shifts on day 13.8 (Fig. [7]) have 
values consistent with other lines in all later observa- 
tions. If these lines ar e shaped by absorpt ion in their red 
wings as proposed by iDrake et al.l ((2008), then the col- 
umn densities have decreased from day 13.8 to day 261. 
The Nvn line on day 26.1 shows an extreme value (also 
in line width, see bottom panel of F ig. [71) , and be l ongs t o 
the new component discovered by INelson et alj (2008); 
however, the velocity measured from the shift of this line 
does not agree with their value of 8, 000-10, 000 km s _1 
derived from the lines between 25- 30 A. While w e mark 
these lines as unidentified in Fig. [31 INelson et alj (2008) 
discuss possible identifications as highly blue-shifted N VI 
and C vi lines. 

The line widths slowly decrease with time. The Nvn 
line at 24.78 A is extremely broad on day 26.1. At this 
time of the evolutio n, this line is part o f the new com- 
ponent reported by INelson et al.1 ((2008?) with a set of 
unidentified lines. There is thus a reasonable chance 
that this line is a blend, making this anomalous veloc- 
ity questionable. Since shock velocities derived from the 
temperatures from the spectral models discussed in 34.21 
represent the evolution of the expansion velocity, they 
can be compared with these values. The shock velocities 
derived from the hottest model component are given in 
the last row of Table [U We discuss the implication of 
the comparison in 



3.2. Line fluxes 

We have used our line fitting program Cora 
(|Ness fc Wichmanf][2002T) to measure line counts which 
are then converted to line fluxes using the effective ar- 
eas extracted from the instrument calibration. The Cora 
program applies the likelihood method described in the 
appendix and adds a model of line templates to the in- 
strumental background. In order to measure the line 
fluxes on top of the continuum, we have added a con- 
stant source background to the instrumental background 
before fitting. The continuum is not expected to change 
significantly over the narrow wavelength range consid- 
ered for line fitting. 

In Table [3J we list the measured line fluxes for the 
strongest lines of Fe and the H-like and He- like lines of N, 
O, Ne, Mg, Si, and S, sorted by wavelength (for emission 
line fluxes measu red on days 39.7, 66.9, and 54 we refer to 
iNess et aLll2007f ). The given uncertainties are the statis- 
tical uncertainties. Additional systematic uncertainties 
arise from the assumed level of the underlying continuum 
and line widths (see H4.5p . The fluxes are not corrected 
for the effects of interstellar or circumstellar absorption. 
In the bottom part of Table [3] we list line ratios of H-likc 
to He-like line fluxes, which are temperature and density 
indicators. 

In the top part of Table [4] we list temperatures de- 
rived from the H-like to He-like line ratios of the same 
species (after correction for Nn as noted in Table Q| us- 
ing theoretical predictions of the same ratios as a func- 
tion of temperature extracted from an atomic database 
computed by the Astrophysical P lasma Emission Code 
(APEC, vl.3: ISmith et al.l l2001al lbh . The temperatures 
are computed under the assumption of collisional equi- 
librium (see end of this section and 34.11 for more de- 
tails) and are average values assuming that the plasma 
is isothermal over the temperature range over which the 
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TABLE 2 

Evolution of line shifts and line widths. 



Ag(A) 




day 13.81 


day 13.88 


day 26.1 


day 39.7 


day 54.0 


day 66.9 


day 111.7 


ID 




HETG 


RGS 


RGS 


LETG 


RGS 


LETGS 


LETGS 


4.73 AA (i 


„A) 


-8.3 ± 1.6 


- 


- 


- 


- 


- 


- 


S XVI a (i 


oA) 


12.7 ± 1.7 


- 


- 


- 


- 


- 


- 


■"shift (km S 


-1j 


-526 ± 103 


- 


- 


- 


- 


- 


- 


'-'width (kms 


-lj 


806 ± 111 


- 


- 


- 


- 


- 


- 


6.18 AA (i 


oA) 


-8.1 ± 0.6 


-4.8 ± 5.7 


- 


-16.1 ± 5.2 


- 


- 


- 


Si xiv <j(i 


uA) 


16.3 ± 0.6 


< 12.1 


- 


17.7 ± 8.9 


- 


- 


- 


""shift (km 8 




-395 ± 29 


-233 ± 275 


- 


-781 ± 253 


- 


- 


- 


''width (kms 


-lj 


790 ± 27 


< 586.2 


- 


861 ± 433 


- 


- 


- 


8.42 AA (i 


nA) 


-21.3 ± 0.7 


-22.4 ± 1.5 


-6.5 ± 2.9 


-25.9 ± 3.0 


-17.5 ± 0.7 


-16.7 ± 8.8 


2.3 ± 18.0 


MgXH <7(i 


oA) 


26.5 ± 0.6 


23.0 ± 3.0 


20.3 ± 5.3 


24.0 ± 5.3 


< 17.1 


< 11.0 


< 27.7 


"shift (km S 


_1 ) 


-759 ± 24 


-796 ± 54 


-230 ± 102 


-924 ± 107 


-623 ± 25 


-593 ± 313 


82 ± 640 


"width (kms 


_1 ) 


945 ± 21 


818 ± 108 


724 ± 189 


853 ± 187 


< 607.8 


< 390.6 


< 985.1 


10.23 AA (i 


a A) 


-14.2 ± 2.3 


- 


- 


- 


-2.2 ± 8.4 


-17.9 ± 10.7 


-1.3 ± 16.4 


Ne X a (r 


nA) 


14.5 ± 2.6 


- 


- 


- 


< 11.3 


< 25.7 


< 12.1 


■"shift (km S 


_1 ) 


-416 ± 69 


- 


- 


- 


-64 ± 246 


-525 ± 314 


-38 ± 482 


"width (kms 


_1 ) 


425 ± 75 


- 


- 


- 


< 331.7 


< 753.7 


< 355.1 


12.13 AA (i 


a A) 


-29.1 ± 1.9 


-17.6 ± 2.4 


-6.1 ± 2.5 


-19.3 ± 2.6 


-3.5 ± 3.3 


- 


-6.8 ± 5.6 


Nex c (r 


oA) 


34.0 ± 1.6 


38.0 ± 3.7 


24.3 ± 4.4 


27.3 ± 4.0 


29.7 ± 5.6 


- 


< 19.0 


"shift (km s 


_1 ) 


-719 ± 47 


-435 ± 60 


-151 ± 62 


-477 ± 65 


-87 ± 81 


- 


-168 ± 139 


"width (kms 


_1 ) 


841 ± 39 


939 ± 91 


601 ± 110 


675 ± 99 


733 ± 137 


- 


< 469.3 


15.01 AA (i 


a A) 


-14.8 ± 3.3 


12.7 ± 4.8 


4.1 ± 3.2 


- 


- 


- 


18.3 ± 4.5 


Fe xvn a (i 


«A) 


19.7 ± 4.4 


65.5 ± 5.8 


30.5 ± 4.4 


- 


- 


- 


28.3 ± 3.0 


"shift (km s 


_1 ) 


-296 ± 65 


254 ± 96 


82 ± 64 


- 


- 


- 


365 ± 91 


"width (km s 


_1 ) 


394 ± 89 


1308 ± 115 


609 ± 89 


- 


- 


- 


565 ± 60 


16.78 AA (i 


qA) 


-19.8 ± 7.3 


-29.6 ± 6.5 


- 


- 


- 


- 


-38.9 ± 8.5 


Fe xvn a (i 


oA) 


15.1 ± 7.3 


13.3 ± 11.2 


- 


- 


- 


- 


< 14.1 


"shift (km s 


_1 ) 


-354 ± 131 


-529 ± 116 


- 


- 


- 


- 


-695 ± 151 


"width (km s 


-1 ) 


270 ± 130 


238 ± 201 


- 


- 


- 


- 


< 252.6 


18.97 AA (i 


qA) 


-79.7 ± 9.6 


-82.6 ± 2.8 


-38.1 ± 5.2 


- 


- 


- 


-26.3 ± 3.0 


O VIII a (r 


oA) 


57.2 ± 5.7 


53.5 ± 3.1 


68.3 ± 3.0 








< 1.0 


"shift (km s 


-1* 


-1260 ± 151 


-1305 ± 44 


-602 ± 82 








-416 ± 47 


"width (km s 


-lj 


904 ± 91 


845 ± 50 


1080 ± 48 








< 16.1 


24.78 AA (i 


qA) 




-106.9 ± 6.3 


30.0 ± 7.7 








-6.9 ± 0.4 


NVII C (r 


nA) 




82.2 ± 5.9 


170.3 ± 1.5 








< 7.9 


"shift (km s 


_1 1 




-1293 ± 76 


363 ± 93 








-83 ± 4 


"width (km s 


-lj 




994 ± 72 


2060 ± 18 








< 95.7 


"shock (Tl) (kms~ 


l,i 


1899 ± 22 


1284 ± 30 


1269 ± 44 






792 ± 37 


"shock(72) (kms- 


1,6 


792 ± 


8 


792 ± 8 


714 ± 7 




747 ± 60 


548 ± 61 



"rest wavelengths 6 From kTi and VT2 in APEC models 




50 100 150 200 250 

day 



Fig. 8. — Evolution of average temperatures derived from ratios 
of H-like and He-like lines for the elements given in the legend. 
The values are listed in Table [J] Error bars are statistical errors; 
systematic errors are small (see £|4.5| |. 

respective H-like and He-like lines are formed. Lines 
originating from high-Z elements probe hotter plasma. 
Since the lines involved in each ratio are from the same 
element, these ratios yield temperature s independently 
of th e elemental abund ances (see, e.g.. ISchmitt fe Nesa 
I2004I ; I Ness et ail |2005) . A graphical representation is 



given in Fig. [5J and it can be seen that the Si ratios 
probe hotter plasma than the N ratios. Significantly dif- 
ferent temperatures are derived, indicating that we are 
dealing with a wide range of temperatures in the early 
observations (before day ~ 70). The measurements for 
the later observations deliver consistent temperatures (at 
least within the large uncertainties), indicating that the 
plasma could be characterized by a single temperature by 
that time. The temperatures derived from the Mg lines 
indicate that the hotter plasma cools until day ~ 70 and 
remains constant after that time. A similar behavior can 
be concluded from the other values but it is not as clear. 
For example the temperatures derived from the Ne lines 
yield a slight increase, however, the temperature on day 
26.1 may also be anomalously low, o wing to line blends 
of the Ncix lines (jNess et al.ll2003al ). We note that the 
N lines observed on day 26.1 show a peculiar behavior in 
the line shifts and -widths (see Tabled, and the derived 
temperature may thus not represent the same plasma as 
the values for the other times of evolution. For days 39.7- 
66.9 the cool component cannot be probed because the 
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O and N lines are outshone by the SSS emission from 
the WD. The temperature derived from the O lines for 
day 3 9.7 is based on the fluxes measured by iNess et al.l 
(|2007f ) on top of the SSS continuum. These fluxes may be 
contaminated by photoexcitations, but the derived tem- 
perature agrees well with the temperatures derived for 
the other days. 

In the bottom part of Table [4] we list densities derived 
from the He-like forbidden-to-intercombination (f/i) line 
ratios for the ions Ovn and Neix. These values have 
been derived assuming collisional equi librium accord- 
ing to the parameterization derived bylGabricl & Jordanl 
(1969), neglecting UV radiation; see, e.g., INess et all 
(2004) for details. The method explores density- 
dependent excitations out of the upper level of the f line 
(ls2p 3 S) into that of the i line (ls2p 3 P). In the low- 
density limit, all ions in the ls2p 3 S state radiate to the 
ground (Is 2 1 S), giving rise to the f line, while with in- 
creasing density, collisional excitations from the Is 2 3 S 
state into the ls2p 3 P state reduce the f line and increase 
the i line. However, the ls2p 3 P-ls2p 3 S transition can 
also be induced by UV radiation, who se presence would 
mimi c a high density if neglected (jBlumenthal et all 
Il972t INess et all 120011) . Especially for the early spec- 
tra we must assume that significant UV contamination 
from the WD has to be accounted for. While the UV 
intensity can b e estimated from I UE observations of the 
1985 outburst (Sh ore et al.lll996l ). the distance between 
the X-ray emitting plasma and the UV source needs to 
be known in order to quantify the contaminating effects 
for the density diagnostics. Since UV radiation fields, if 
present, mimic high densities, we treat the values with 
great caution, but we can at least conclude that the den- 
sity is not higher than any of the values listed in Tabled) 

In the next section we present models for which the 
assumption of optically thin plasma is made. As a test 
of this assumption we measure the line flux ratio of the 
Ly/3 (3p-ls) to Lya (2p-ls) lines of the H-like ion Nex 
(10.23A and 12.12A, respectively), which increases 
with increasing optical depth, but also depends on the 
electron temperature and the amount of photoelectric 
absorption in the line of sight. Theoretical predictions of 
the same ratio for an optically thin plasma in collisional 
equilibri um, from the ato mic databases CHIANTI ver- 
sion 5.2 (|Landi et al.ll2006D and APEC, vary between 0.1 
and 0.3 within the temperature range 10 6 to 10 s K, while 
for day 13.8 we measure a ratio of 0.17 ± 0.02. This is 
well within th e expected ran g e. We also refer to the tests 
presented by I Nelson et al.1 (|2008f ) who measured the 
so-called G-ratio of the He-like triplets of S, Si, and Mg 
and concluded that the plasma is collisionally dominated. 



4. ANALYSIS 

For the interpretation of the observations described 
above we use two model approaches. First we com- 
pute multi-tem perature plasm a models with the fitting 
package xspec (jArnau d 199G). We use the atomic data 
computed by APEC which are similar to the MEKAL 
database, and the r e sults can be compared to those 
given by Bo de et"aT1 (|2006h for earlier X-ray observa- 
tions. Next, we use the measured emission line fluxes 
from Table to construct a model of a smooth temper- 
ature distribution. This model allows us to determine 



relative abundances. 

4.1. Description of methods and model assumptions 

For our modeling of the X-ray spectra of the shock we 
assume: (1) that all emission originates from the same 
volume with the same abundances, (2) that the plasma 
is in a collisional equilibrium, and (3) that it is optically 
thin. 

The first assumption is implicit in all spectral analyses 
in X-rays unless spatial resolution is available. Although 
there is likely stratification to some extent in the emit- 
ting environment, we have no basis on which we can de- 
velop more refined models. We further have to assume a 
uniform plasma, which implies that interstellar and cir- 
cumstellar absorption can only be modelled with a single 
absorption component. 

The line ratios used in the end of §3.21 support the 
other two assumptions. In a collisional plasma, all tem- 
peratures are kinetic temperatures, derived from the dis- 
tribution of velocities, which is commonly assumed to be 
Maxwellian. While in a shocked plasma collisions are 
the main energy source for all atomic transitions, the as- 
sumption of an equilibrium is not necessarily valid. Note 
that rapid recombination can lead to non-equilibrium 
conditions, since recombination into excited states leads 
to an overpopulation of upper levels and consequently to 
excessively high fluxes in certain emission lines. Never- 
theless, we base our analysis on the assumption of equi- 
librium conditions and discuss the implications of this 
assumption where relevant. 

The third assumption implies that all emission that 
is produced in the collisional plasma escapes unaltered. 
However, lines with high oscillator strengths may be re- 
absorbed within the plasma and reemitted in a differ- 
ent direction (resonant line scattering). Depending on 
the plasma geometry, resonance lines may be stronger 
or weaker compared to optically thin plasma, since res- 
onance line photons can be scattered out of the line of 
sight or into the line of sight. In a spherical geometry 
the processes cancel out and no effects from resonant 
line scattering are detectable. Ways to detect res onant 
line scattering are discussed bv lNess et ah ( 2003b). The 
emission line fluxes presented in §3.21 show no signatures 
of resonant line scattering. 

In a collisional plasma the brightness of a source is ex- 
pressed in terms of the volume emission measure, VEM, 
which is a measure of the intensity per unit volume (in 
cm~ 3 ). The VEM is defined as VEM = J n 2 c dV with 
n e the electron density and V the emitting volume. A 
given value of VEM is thus proportional to the emitting 
volume; however, volumes can only be determined from 
independent density measurements which are difficult to 
obtain from X-ray spectra (see M3.2|) . 

The volume emission measure as a function of tem- 
perature, T, is called the emission measure distribution 
(EMD). A given EMD is a model that allows the calcula- 
tion of continuum emission by bremsstrahlung and emis- 
sion line fluxes, which together form a predicted X-ray 
spectrum that can be compared to an observed spectrum. 
In order to predict the line flux for a line at wavelength A, 
one needs the respective line contribution function with 
temperature T 



A ni(T)n e (T) he 47r<i 2 
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TABLE 3 
Evolution of emission line fluxes. 







flux 6 


flux 6 


flux 6 


flux 6 


flux 6 


flux 6 


flux 6 


Ion 


(A) 


day 13.81 


day 13.88 


day 26.1 


day 111.7 


day 205.9 


days 204-208 6 


day 239.2 


Fe xxvi 


1.78 


< 144 








_ 


_ 




Fe xxv c 


1.85 


1037 ± 190 








_ 


_ 


_ 


Sxvi 


4.73 


219 ± 22 




day 39.7: < 33.8 


— 


_ 




Sxv c 


5.04 


624 ± 84 




day 39.7: 


78.1 ± 67.7 


_ 


_ 




Si xiv 


6.18 


349 ± 13 


260 ± 32 


138 ± 29 


2.7 ± 1.8 




< 1.5 




Si xiii c 


6.65 


620 ± 32 


795 ± 112 


560 ± 124 


10.3 ± 5.6 


_ 


< 5.1 


_ 


Mgxn 


8.42 


412 ± 11 


585 ± 19 


269 ± 19 


5.1 ± 2.0 


< 0.3 


1.3 ± 0.7 


< 0.3 


Mg xi c 


9.20 


383 ± 30 


369 ± 35 


374 ± 47 


18.8 ± 7.6 


4.5 ± 2.9 


4.1 ± 2.4 


2.2 ± 1.1 


Nex 


12.13 


273 ± 20 


309 ± 9.2 


309 ± 12 


21.6 ± 3.0 


3.0 ± 0.7 


3.2 ± 0.7 


2.0 ± 0.5 


+Fe xvn 


12.12 


1.08 times the flux at 12.26 A if logT < 6. 


9 otherwise ne; 


jligiblc 






Fe xvn 


12.26 


52.8 ± 6.8 


52.5 ± 6.6 


82.8 ± 9.3 


3.2 ± 1.8 


< 1.6 


1.5 ± 0.6 


< 0.75 


+Fe xxi 


12.28 


at logT > 6 


.9 only Fe xxi otherwise only Fe XVII 








Ne ix 


13.44 


71.8 ± 8.9 


99.8 ± 5.3 


148 ± 9.4 


19.6 ± 2.8 


3.7 ± 0.8 


3.4 ± 0.7 


3.1 ± 0.5 


Neix (i) d 


13.55 


19.8 ± 6.8 


38.2 ± 4.5 


60.3 ± 8.4 


3.4 ± 1.8 


1.2 ± 0.6 


1.0 ± 0.6 


0.7 ± 0.4 


Ncix (f) d 


13.69 


33.9 ± 6.8 


53.3 ± 4.1 


82.3 ± 7.4 


10.1 ± 2.1 


3.3 ± 0.7 


2.8 ± 0.6 


2.4 ± 0.5 


Fe xvn 


14.21 


33.6 ± 7.5 


45.5 ± 3.4 


59.5 ± 5.6 


1.8 ± 1.4 


1.6 ± 0.6 


1.0 ± 0.5 


- 


Fe xvn 


15.01 


67.3 ± 8.1 


69.5 ± 3.6 


156 ± 7.7 


14.3 ± 2.1 


3.4 ± 0.7 


2.7 ± 0.5 


2.7 ± 0.5 


O vni 


18.97 


84.2 ± 19 


92.5 ± 3.1 


186 ± 13 


33.5 ± 2.8 


7.7 ± 1.4 


5.3 ± 0.6 


5.2 ± 0.5 


Ovn 


21.60 




14.8 ± 1.9 


29.5 ± 3.1 


6.8 ± 1.8 


2.4 ± 0.5 


2.1 ± 0.6 


0.9 ± 0.3 


Ovu (i) d 


21.80 




< 1.5 


13.6 ± 2.4 


1.3 ± 1.3 


1.4 ± 0.4 


< 0.5 


0.6 ± 0.3 


Ovn (f) d 


22.10 




6.4 ± 1.3 


33.4 ± 3.2 


5.4 ± 1.7 


2.0 ± 0.5 


0.8 ± 0.6 


1.0 ± 0.3 


Nvn 


24.78 




43.0 ± 1.8 


259 ± 12.8 


15.7 ± 2.3 


6.6 ± 0.7 


3.9 ± 0.6 


4.0 ± 0.5 


Nvi 


28.78 




2.3 ± 0.6 


35.4 ± 3.5 


1.8 ± 1.1 


< 0.4 


< 1.5 


0.4 ± 0.2 


Cvi 


33.74 




< 0.6 


< 0.5 


< 1.4 


< 0.3 


< 0.9 


< 0.1 


Temperature-sensitive line ratios 














Fexxvi/Fexxv 




< 0.2 








- 


- 


- 


Sxvi/Sxv 




0.35 ± 0.06 




day 39. 


7: < 3.25 








Sixiv/Sixm 




0.56 ± 0.04 


0.33 ± 0.06 


0.25 ± 0.08 


0.26 ± 0.23 








Mgxn/Mgxi 




1.08 ± 0.09 


1.59 ± 0.16 


0.72 ± 0.10 


0.27 ± 0.15 


< 0.2 


0.32 ± 0.25 


< 0.3 


Neix/Nex 




3.80 ± 0.54 


3.09 ± 0.19 


2.10 ± 0.16 


1.10 ± 0.22 


0.81 ± 0.26 


0.94 ± 0.28 


0.65 ± 0.19 


Oviii/O vii 






6.25 ± 0.83 


6.31 ± 0.79 


4.93 ± 1.37 


3.21 ± 0.89 


2.52 ± 0.78 


5.78 ± 2.00 


N vo/N vi 






18.70 ± 4.94 


7.32 ± 0.81 


8.72 ± 5.48 


> 14.8 


> 2.2 


10.00 ± 5.15 


Density-sensitive line ratios 














Neix (f/i) 




1.7 ± 0.7 


1.4 ± 0.2 


1.4 ± 0.2 


3.0 ± 1.7 


2.8 ± 1.5 


2.8 ± 1.8 


3.4 ± 2.1 


O vn (f/i) 






> 4.3 


2.5 ± 0.5 


> 4.2 


1.4 ± 0.5 


> 4.0 


1.7 ± 1.0 



Uncertainties and upper limits are 68.3%. Additional systematic uncertainties are discussed in £14.51 • "rest wavelengths 

• 6 10 -14 ergcm" 2 s" 1 • c sum of three lines • d intersystem line • e Sum of Chandra spectra taken between days 203.3 and 208.2 (see 

Tabled}. 



TABLE 4 

Temperatures and Densities derived from H-like to He-like line ratios 

Ion day 13.81 day 26.1 day 39.7 day 54 day 66.9 day 111.7 days 204-208° day 239.2 

logT in K, after correction of line ratios for ./Vh with values 5 X 10 21 cm - 2 for days 13.8 and 26.1 and 2.4 X 10 21 cm -2 for the rest. 
Fe < 7.8 

S 7.2lt°'° 2 < 7.61 

Si 7.09 ± 0.01 6.98 ± 0.03 7.02 ± 0.02 7.01+g ° 5 7.05±q °1 7.00±g' °9 

Mg 6.96 ± 0.01 6.91 ± 0.02 6.87 ± 0.02 6.85 ± 0.02 6.81 ± 0.03 6.80±°;° 4 6.8ljl° % < 6.8 

Ne 6.81 ± 0.02 6.72 ± 0.01 6.74 ± 0.02 6.77 ± 0.03 6.79i°;f 6.65 ± 0.03 6.63 ± 0.04 6.59 ± 0.03 

O 6.46 ± 0.02 6.46 ± 0.02 6.42 ± 0.03 6.52 ± 0.04 6.43 ± 0.04 6.55±o° fi 

N 6.36 ± 0.04 6.23 ± 0.01 - 6.40l u ;°g > 6.5 6.42 ± 0.11 

Densities, logra,. in cm~ 3 . from He-like triplet ratios assuming no UV illumination (see text £13-2 1) 

Ovn <8 10.3l ( ( j;2 < 8 10 - 7 -o% 10 - 6 ±o:s 

Neix 11.7 ± 0.4 H-9t°;4 ~ > 10.3 > 10.8 >9 
"Sum of Chandra spectra taken between days 203.3 and 208.2. The results arc consistent with XMM observations taken day 205.9. 

with h Planck's constant, c the speed of light and A the tion, and n e is the electron density. The ratio ni(T)/riE 

wavelength of the line. The number densities are given is the ionization balance, and calculations of nj(T)/riE 

with n with subscripts u for ions in the upper level in the as a function of tem perature have been presented by 

transition, / for the ionization stage in which the transi- iMazzotta et al.l (|1998f) . The ratio tie /nn represents the 

tion occurs, E for the element giving rise to the transi- absolute elemental abundance, and A u i is the Einstein 
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A-value. With a given EMD (which can also be written 
as V EM{T)) the predicted line fluxes are 

F\ = J VEM(T) x G\(T) x dT. (2) 

An EMD can be defined as one or more isothermal com- 
ponents ( fl4.2p or as a continuous function of temperature 

When fitting models using XSPEC f ^4.2[) all informa- 
tion that is available in the atomic database is used to 
constrain the models. We make use of a second approach 
f q4.3[) selecting only the most reliable atomic data, and 
constrain the models by fitting the measured line fluxes 
rather than the entire spectrum. Each approach has 
strengths and weaknesses, and we compare the two ap- 
proaches in i )4.41 

In §4.5l wc discuss estimates of systematic uncertainties 
in addition to the given statistical uncertainties. 

4.2. Multi-temperature plasma models 

We construct spectral models u sing the X-ray fitting 
package xspec version 11.3.2ag (|Arnaudl 119961 ) which 
combines all information in a given atomic database and 
generates count spectra to be compared to the observed 
count spectra. We chose multi-temperature APEC mod- 
els comprising the sum of n independent isothermal com- 
ponents with v ariable abundances relative to the solar 
values listed by lAnders fe; Grevessd (|1989f ). Line fluxes 
are obta ined applying Eq. [2] u sing the ionization bal- 
ance by iMazzotta et al.l (1998), and the plasma den- 
sity is assumed constant at logn e = (jSmith et al.1 
l2001al) . We vary only the abundances of elements which 
produce strong emission lines in the respective spec- 
tra (see Figs. [2] to [5]) and assume solar abundances for 
the other elements. We correct for interst ellar absorp- 
tion u sing the tbabs module developed bv lWilms et al.l 
(2000) and allow the neutral hydrogen column density, 
Ah, to vary only for the first two observations. For the 
later observations we fix Ah at the interstellar value of 
Ah = 2.4 x 10 21 cm~ 2 . The effects of Ah have a stronger 
influence on the cooler component and may affect the 
abundances of elements whose lines are formed at longer 
wavelength (i.e., nitrogen and oxygen). 

We use the optimization procedures provided by the 
xspec programme to obtain best fits and the XSPEC 
command error to calculate the parameter uncertain- 
ties, yielding 90-per cent uncertainties. This command 
steps through a range for a given parameter, and in the 
course of this process further improvements of the fit 
can be found. The uncertainties returned by the error 
command are only statistical uncertainties that describe 
the precision of the measurement but not necessarily the 
accuracy of the respective parameters (see fc|4.5|) . The 
results are listed in Table [3 and some of the correspond- 
ing models are shown in Fig. [9] The elemental abun- 
dances are only varied for day 13.8, because this is the 
best d ataset. The values relative to lAnders &: Grevessel 
(1989) adopted for all observations are given in Table El 
middle column. 

While the true EMD is most likely a continuous distri- 
bution, we use multi-temperature models. The only con- 
tinuous EMD models to chose from in XSPEC are Cheby- 
shev polynomials of no more than six orders and are not 
constrained to be positive. We regard these EMDs as 



not sufficient for our purposes and therefore use the more 
standard multi-T models. The number of free parame- 
ters, and thus the number of temperature components, 
has to be chosen to be as small as possible while still 
achieving a good fit. The spectra taken on day 13.8 
have high statistical quality, and a 3-temperature (=3-T) 
model yields significantly better reproduction of the data 
than 2-temperature models. With variable abundances, 
no fourth temperature component is required to improve 
the fit. The APEC model has a redshift parameter that 
we allow to vary in order to account for line shifts (see 
§3.ip . however, we cannot account for line broadening 
in excess of the instrumental line broadening function. 
While the model could be folded with a Gaussian with 
variable width, this is computationally expensive and un- 
feasible with the given high number of spectral bins and 
free parameters to be iterated. The long- wavelength lines 
and the associated elemental abundances (particularly of 
N and O) may thus be poorly determined, and our sec- 
ond approach ( £ 14. 3[) is more reliable for the abundances 
of these elements. Meanwhile, the lines at shorter wave- 
lengths are not broadened by as much, and the abun- 
dances of the other elements are less affected. We fit the 
model simultaneously to the HEG and MEG spectra plus 
the RGS1 and RGS2 spectra (top panels in Fig. [9]). Our 
model agrees bett er with the RGS dat a (top panel) than 
that presented bv iNelson et al.1 (|2008| ). fig. 3, who have 
not varied the abundances but used four temperature 
components. No formal value of \ 2 was given, but visual 
inspection clearly shows that their model does not re- 
produce the RGS spectrum. This demonstrates that the 
effects from non-solar abundances are indeed detectable. 

Since on day 26.1 the source was still bright, the ob- 
served spectra are also of high statistical quality, and 3-T 
models are better than 2-temperature models. Since we 
expect no detectable changes in the composition, we use 
the elemental abundances found from the day 13.8 obser- 
vations for this and the following datasets. We discard all 
spectral bins longward of 20 A because the emission does 
not originate from the shock (third and fourth panels in 
Fig. EI) . The N abundance is now less certain because the 
strong Nvn line at 24.78 A is excluded. We concentrate 
on the RGS spectra but also compute a model including 
the MOS1 spectra which are sensitive at higher energies 
and are thus suited to constrain the hotter component. 
As can be seen from Table [5J the model parameters are 
identical, and only the uncertainties of the model with 
the MOS1 data included are smaller. The hot component 
is thus detectable with the RGS alone. 

Since the spectra taken on days 39.7-66.9 are com- 
promised by the SSS emission (bottom panel of Fig. [9]) 
the X-ray emission from the low-temperature shocked 
plasma cannot be probed. The spectrum shortwards 
of ~ 14 A is not well enough exposed to require three 
temperature components. However, in order to compare 
the results, we allow three temperature components with 
the option that the fitting procedure can assign small 
values of emission measure to those components that are 
not detectable. We caution, however, that the hottest 
temperature component that probes the bremsstrahlung 
continuum may be overestimated due to systematic 
uncertainties in the instrumental background (see tj4.5[) . 
If the theoretical continuum is of the same order as the 
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TABLE 5 

Model parameters from multi-T APEC models 





day 13.8 


day 26. 1 6 


day 26.1 c 


day 39.7 


day 66.9 


day 111.7 


day 239.2 




4.19 ± 0.11 


1.80 ± 0.24 


1.91 ± 0.07 


1.88 ± 0.13 


74.4l| 2 5 


0.74 ± 0.06 


< 79.9 


log(T)° 


7.69 ± 0.01 


7.32 ± 0.06 


7.35 ± 0.02 


7.34 ± 0.03 


8.931°;?* 


6.93 ± 0.04 


< 8.97 


log(EMi) a 


58.02 ± 0.01 


57.48 ± 0.04 


57.51 ± 0.02 


57.23 ± 0.03 


56.6 ± 0.14 


56.15 ± 0.09 


< 54.04 




0.74 ± 0.01 


0.74 ± 0.02 


0.74 ± 0.01 


0.60 ± 0.01 


66+ 13 

u - DD -0.09 


0.35 ± 0.03 


0.37 ± 0.02 


log(T) a 


6.93 ± 0.01 


6.93 ± 0.02 


6.93 ± 0.01 


6.84 ± 0.01 


fi SS+0-08 


6.61 ± 0.11 


6.63 ± 0.03 


log(EM 2 ) a 


57.76 ± 0.01 


57.69 ± 0.05 


57.67 ± 0.04 


57.44 ± 0.02 


56.81 ± 0.20 


55.21 ± 0.09 


55.65 ± 0.03 


kT 3 


0.30 ± 0.01 


0.37 ± 0.03 


0.38 ± 0.02 


0.11 ± 0.01 


0.14 ± 0.06 


0.12 ± 0.03 


0.11 ± 0.02 


log(T) a 


6.55 ± 0.01 


6.64 ± 0.96 


6.64 ± 0.03 


6.10 ± 0.06 


6.21 ± 0.24 


6.14 ± 0.11 


6.12 ± 0.10 


\og(EM 3 ) a 


57.46 ± 0.02 


57.60 ± 0.09 


57.61 ± 0.05 




58.88 ± 1.21 


55.88 ± 0.25 


55.32 ± 0.45 


N H 


6.95 ± 0.35 


5.59 ± 0.02 


5.56 ± 0.13 


2.4 


2.4 


2.4 


2.4 




0.67, 21317 


1.70, 2939 


1.73, 3151 


(0.65, 1928) e 


(0.18, 1294) e 


(0.21, 10088) e 


(0.59, 5228) e 



"Units: kT in koV, T in K, EM in cm" 3 , N a in 10 21 cm" 2 • 6 Only RGS data • C RGS and MOS1 data simultaneously • ^degrees of 
freedom • e after iteration with C-statistics • * fixed at values from day 13.8 



noise in the background, arbitrarily high temperatures 
may result which will have to be treated with caution. 
Since the spectra taken after day 26.1 c ontain man y 
bins with low counts, we use C-statistics (|Cashl [l979h . 
This approach is the based on the maximum likelihood 
(ML) method described in the appendix. We calculate 
a formal value of x^ od after fitting with est at for 
comparison with the other fits. We use the errors on the 
count rates from the extracted spectra. 

Because the spectra on days 13.8 and 26.2 are of such 
high quality, we investigate absolute abundances using 
these two datasets (Fig.[TU]). Although no hydrogen lines 
are present in the X-ray range, the absolute abundances 
can be determined from the strength of the continuum 
relative to the lines. The brightness of the continuum de- 
pends on the number of free electrons which, in an ionized 
plasma, scales with the hydrogen abundance. We thus 
need spectra with sufficient continuum emission. We step 
through a grid of (fixed) abundances and fit the remain- 
ing parameters to minimize \ 2 ■ I n Fig- EHl we show the 
relative changes in y 2 for each grid point in compari- 
son with the 68-% and 95-% confidence ranges (for 14 
free parameters). While for day 13.8, solar abundances 
are preferred, the spectrum taken on day 26.1 suggests 
a somewhat lower metallicity, but from the confidence 
intervals one can see that their determination is highly 
uncertain. We therefore fix the absolute abundance at 
solar values and concentrate on the relative abundances. 

The final model parameters are summarized in TableO 
From top to bottom we list log T and log VEM for each 
component, value of TVh, and Xred with number of degrees 
of fre edom (dof). The element al abundances relative to 
solar ([Anders fc Greves sc 1989) as determined from the 
day 13.8 dataset are given in Table [7] and have been used 
for all models. 

In the same way as iBode et all (|2006l ) we compute 
shock velocities, v s hock, from the temperatures of the 
first and second component of the APEC models (see 
Table [S]), and from the highest temperatures found from 
line ratios (Table 2]). In Fig. [TT1 we com pare these re- 
sults with those from IBode et all (|2006f ). The dotted 
and dashed lines indicate the expected evolution for a 
radiatively cooling plasma and an adiabatic plasma, re- 
spectively (jBode fc Kahnlfl985T) . The velocities derived 
from Ti of our 3-T APEC models for days 13.8 and 26.1 



are slightly higher than the Bo de values. The re ason 
is that the 1-T models used by IBode et all (|2006l ) are 
an average of all temperature components, accounting 
for some of the cooler plasma that in our 3-T models 
are accounted for by the two cooler components. The 
evolution o f T\ fo llows the same trend as observed by 
IBode et al.l (|2006f ). After day 26.2, the hottest compo- 
nent is much fainter, and T\ is less certain. In Figs. [4] 
and [5] one can see that the continuum emission level is 
significantly lower than that seen in Figs. [5] and [31 Since 
the parameters of the hottest temperature component 
arc dominated by the continuum, systematic uncertain- 
ties from background noise have a stronger effect, and 
the velocities derived from the hottest temperature com- 
ponents may be overestimated for the observations taken 
after day 26.1 (see gTS]) . 

The second plasma component is very similar to the 
values derived from the line ratios (Tabled]), but these 
curves follow a different trend than the hottest compo- 
nent. For the observations of days 13.8 and 26.1, the 
line ratios yield much lower velocities than those from 
the APEC models. We attribute this difference to the 
stronger continuum observed for these two days. Since 
the continuum is dominated by the hottest plasma, the 
hottest temperature in the APEC models is driven by 
the continuum which, during the early observations, re- 
flects a higher temperature than any of the emission lines 
can probe. Meanwhile, the emission lines can probe the 
structure of the temperature distribution better, and in 
the next section we describe an approach that focuses on 
a few selected emission lines. 

4.3. Emission measure modeling 

We use the measured line fluxes listed in Table [3] (cor- 
rected for absorption; see below) in order to reconstruct a 
continuous mean emission measure distribution (EMD) 
as a function of temperature, i.e., VEM(T). We as- 
sume a constant electron pressure of log(P c ) = 13.0 
(in units Kcm~ 3 ), which is equivalent to a density 
log(n c [cm~ 3 ]) < 7, depending on temperature. While 
this assumption is more realistic than log(n e ) = (as 
used for the APEC models) this is still very crude, how- 
ever, we use only lines that are not density-sensitive such 
that the results do not depend on the assumed pressure. 

We concentrate on the two simultaneous observations 
taken on day 13.8 because the combined Chandra and 
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Fig. 9.— Best-fit APEC models (from top to bottom) to the XMM-Newton/RGS and Chandra/MEG spectra on day 13.8 XMM- 
Newton /MOS1 and RGS spectra taken on day 26.1, and the Cfeancira/LETG spectrum taken on day 39.7. For day 26.1 the soft component 
longwards of 20 A was excluded from the fit, and for day 39.7 the SSS spectrum had to be excluded. 

120 C ' ' ' ' 
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0.4 0.6 0.8 

abundance (rel. to solar) 

Fig. 10. — Elem ental abundances relative to solar 

(Anders & Grcvcssc 1989) from APEC fits to the spectra of days 
13.8 and 26.1. Each grid point is the result of iterating all free pa- 
rameters with the absolute abundance fixed at the respective grid 
value. Plotted are the changes in x 2 compared to the best fit. The 
dashed horizontal lines mark the 68-% and 95-% confidence ranges. 

XMM-Newton data provide the largest coverage in lines 
and most reliable line fl ux measurement s . Det ails of the 
method are described in iNess &: Jordan! (|2008h . A simi- 



lar approach has been described by INess et alj (|2005f) . 
A continuous EMD is a more realistic representation 
than multi-temperature models, and it is a way to over- 
come the extreme simplification of a ssuming an isother- 
mal plasma (e.g., iBode et aD I2006T ) . We use Eq. [2] to 
compute line fluxes from a given EMD and compare the 
predicted fluxes to the measured fluxes. For Eq. [J we 
assum e the same ionization balance that INess &: Jordan! 
(2008) used. The elementa l abun dances are relative to 
solar bv lGrevesse fc Sauvall (| 19981 ). 

As a guide to construct a starting EMD we compute 
the so-called emission measure loci L\ (T) , which are the 
ratios of the measured line fluxes, fx, and the line con- 
tribution functions, G\(T) (Eq. [T|), i.e., 

^> = awy (3) 

In Fig. [12] we show these loci for a set of lines selected 
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Fig. 11. — Evolutio n of the shock vel ocity obtained from the 
temperatures found by Bode ct al. (2006), from the temperatures 
of the first and second components of the APEC models (Table O, 
and the highest temperatures measured from line ratios (Table [1J. 
The dashed and dotted lines indicate the expected power law decay 
in an adiabatic plasma (f 1 ^ 3 ) and in a radiatively cooling plasma 
(t~ 1 ' 2 ), respectively. The erro r ba rs are statistical errors. For 
additional systematic errors see £|4,5I 



TABLE 6 

Comparison of line flux measurements with predictions 

FOR DAY 13.8 



Ion" 


A 


F b 

mcas 


fprcd 


N VI 


28.79 


543 


0.93 


± 


0.23 


N VII 


24.78 


1514 


1.03 


± 


0.04 


O VII 


21.60 


1200 


1.13 


± 


0.14 


ovii 


22.10 


686 


1.00 


± 


0.20 


O VIII 


18.97 


1871 


0.86 


± 


0.19 


Ne IX 


13.45 


358 


0.95 


± 


0.12 


NelX 


13.70 


183 


0.96 


± 


0.19 


NeX 


12.14 


927 


0.99 


± 


0.07 


Mg XI 


9.17 


705 


1.03 


± 


0.08 


Mg XII 


8.42 


669 


0.99 


± 


0.03 


Si XIII 


6.65 


816 


1.05 


± 


0.05 


Si XIV 


6.19 


438 


0.94 


± 


0.03 


SXV 


5.04 


708 


1.10 


± 


0.15 


S XVI 


4.73 


245 


0.72 


± 


0.07 


Fe XVII 


12.26 


65.0 


1.13 


± 


0.61 


Fe XVII 


15.02 


448 


0.85 


± 


0.10 


Fc XVIII 


14.21 


216 


1.17 


± 


0.26 


FcXX 


12.83 


69.8 


0.69 


± 


0.39 


FeXX 


12.85 


97.7 


1.08 


± 


0.44 


FeXXI 


12.28 


174 


0.96 


± 


0.39 


Fe XXII 


11.77 


92.5 


1.07 


± 


0.20 


Fc XXIII 


10.98 


69.6 


1.33 


± 


0.26 


Fe XXIII 


11.74 


99.8 


1.49 


± 


0.25 


Fe XXIV 


11.17 


46.8 


1.46 


± 


0.41 


FcXXV 


1.85 


1045 


29.4 


± 


5.40 


Fe XXVI 


1.78 


< 145 


< 


5.19 



• "Lines used in deriving the mean EMD are given in bold face. 

• ^Fluxes from Table[3]in 10 — 14 erg cm -2 s _1 , corrected for 
absorption assuming Nn = 5 X 10 21 cm -2 • c Predicted from the 
derived EMD, assuming constant pressure logP c = 13.0 



on the grounds that the atomic physics are reliable and 
that the measured line fluxes (corrected for Nh = 5 x 
10 21 cm" 2 for day 13.8 and N H = 2.4 x 10 21 cm" 2 for day 
111.7) are either not blended with other lines or are easy 
to deblend (see comments in Table [3]) . A few other lines 
are shown for comparison in light gray with the label 
at their minima. Since the line contribution functions 
G\(T) scale with the elemental abundances (see Eq. [1]), 
a reduction in abundances leads to an increase of L\ (T) 



TABLE 7 

Elemental abundances relative to solar 





day 13.8 


day 111.7 




EMD model 


APEC model 


EMD model 


correction factors 






O 


1.0 


0.77 ± 0.03 


1.0 


N 


6.94 ± 0.44 


7.12 ± 0.70 


7.54I1 2 


Ne 


1.39 ± 0.11 


0.68 ± 0.03 


1.36 ± 0.16 


Mg 


2.29 ± 0.08 


1.06 ± 0.03 


1.91 ± 0.70 


Si 


2.04 ± 0.08 


0.88 ± 0.03 


1.62 ± 1.0 


S 


3.16 ± 0.21 


0.69 ± 0.12 


— 


Fe 


0.50 ± 0.04 


0.29 ± 0.01 


0.43 ± 0.10 


relative to 


oxygen 






[N/O] 


0.84 ± 0.03 


0.95 ± 0.05 


0.882^q'47 


[Ne/O] 


0.14 ± 0.04 


0.04 ± 0.03 


0.13 ± 0.05 


[Mg/O] 


0.36 ± 0.02 


0.24 ± 0.02 


0.28 ± 0.20 


[Si/O] 


0.31 ± 0.02 


0.16 ± 0.02 


0.21 ± 0.42 


[S/O] 


0.50 ± 0.03 


0.17 ± 0.08 




[Fe/O] 


-0.30 ± 0.03 


-0.50 ± 0.03 


-0.37 ± 0.11 


relative to 


iron 






[N/Fe] 


1.14 ± 0.05 


1.43 ± 0.07 


1.24 ± 0.58 


[O/Fe] 


0.30 ± 0.05 


0.50 ± 0.03 


0.37 ± 0.09 


[Ne/Fe] 


0.45 ± 0.05 


0.53 ± 0.05 


0.47 ± 0.15 


[Mg/Fe] 


0.66 ± 0.04 


0.73 ± 0.03 


065 ± 0.23 


[Si/Fe] 


0.61 ± 0.04 


0.65 ± 0.04 


0.58 ± 0.35 


[S/Fc] 


0.80 ± 0.05 


0.67 ± 0.12 





all values relative to iGr evcssc & Sauva] 1)19981 ) 



at all temperatures and vice versa. Because most Fe lines 
are quite weak, difficult to measure, and are subject to 
less certain atomic physics because of the complex ion 
structures, we exclude all Fe lines from constraining the 
mean EMD. 

To find a model that reproduces the selected line fluxes 
we construct an initial EMD by eye. We start with the 
envelope curve below the minima of all L\ (T) curves and 
make successive changes to the EMD and to the elemen- 
tal abundances until the predicted line fluxes agree qual- 
itatively with the measured values. We adjust only the 
abundances of elements that produce strong lines, ex- 
cept for oxygen. The oxygen line fluxes pose a constraint 
on the normalization, and all abundances are thus rela- 
tive to oxygen. Since the line contribution functions are 
broader than 0.2 dex, no narrow features in the temper- 
ature distribution can uniquely be resolved, and we thus 
allow no features in the EMD that are narrower than 
the line contribution functions. We stress that with our 
approach we are determining only one possible represen- 
tation of the true nature of the shocked plasma since 
Eq. [2] represents a Fredholm inte gral equation which is 
not uniquely solvable. However, iNess fc Jordan! (|2008l ) 
pointed out that the determination of elemental abun- 
dances seems fairly robust against the precise form of 
the assumed mean emission measure distribution. 

Once a reasonable model is found we fine-tune the 
model by iteration of the mean EMD, optimizing the 
predicted line fluxes using the method described by 
INess fc Jordan! (|2008f ). Based on the ratios of measured 
to predicted line fluxes for the best-fit models, we modify 
the abundances of elements where systematic discrepan- 
cies can be identified and repeat the fine-tuning of the 
EMD. In this way we consecutively approach a good rep- 
resentation of all lines included in the fit as well as some 
other lines that are not included. 
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Fig. 12. — Upper panel: Emission measure loci L \(T), calculated for each line using Eq. [3] with line fluxes measured on day 13.8, and 
line emissivities assuming rescaled solar abundances (Grcvcssc & Sauval 1998) using the scaling factors listed in the top part of Tabic 171 
The black bullets, connected by a thick solid line indicate the mean emission measure distribution (EMD) yielding the best reproduction of 
the measured line fluxes. The loci for some lines that are not use d to optimize the EM D are shown with light gray. The black X symbols 
indicate the results from 3-T APEC models listed in Table \5\ (see £|4.2l and discussion in i|4,4l l. Middle panel: Emission measure loci only 
for Fe lines (gray: assuming solar abundances, black: corrected by factor 0.46 times solar) and the best-fit EMD from the top panel in 
purple. The bulk of the lines demands a reduction of the Fe abundance, but the Fexxv locus (top right curve) is then too high. The locus 
of Fexxvi (far right curve) is an upper limit. The contribution function of H-like Fexxvi above 10 s K is estimated by extrapolating from 
the data below 10 s K assuming the same shape as the H-like line of Sxvi. Bottom panel: Same as top panel for line fluxes measured on 
day 111.7. The same elemental abundances are used. The black X symbols indicate the results from APEC models listed in Table[5] 



The final model is indicated with the thick solid line in 
the top panel of Fig. [T2J The best fit yields two peaks at 
~ 10 7 K and ~2x 10 6 K. The high-temperature regime 
is poorly determined because the only lines formed at 
temperatures above logT = 7.5 are those of Fexxv and 
Fexxvi. For the latter line we only have an upper limit 
to the flux. The use of these lines is limited by the un- 
known Fe abundance at this stage. 

In the middle panel of Fig. fTS] we show emission mea- 
sure loci of ten Fe lines measured from the MEG and 
HEG spectra taken on day 13.8. The grey curves are 
the loci calculated with solar Fe abundance, and all loci 
around logT ~ 7 need to be raised (yielding a reduc- 
tion of the Fe abundance) in order to be consistent with 
the mean EMD derived from the other lines. If the Fe 
abundance is reduced by a factor 0.5 (black loci), the re- 



production of the Fe lines improves significantly. Only, 
the Fexxv line is not reproduced, yielding an underpre- 
diction by a factor of 30. 

The excessively high Fexxv flux in combination with 
the non-detection of the Fexxvi line is difficult to ex- 
plain. While the underestimated flux for Fexxv could 
be fixed with more emission measure at high tempera- 
tures, such a modification demands a detectable flux of 
the Fexxvi line. When increasing the emission mea- 
sure only at temperatures where the Fexxv lines are 
formed, the S xvi line is significantly overpredicted. We 
are confident that the Fexxv emissivity function is not 
underestimated, as the two atomic data bases APEC and 
CHIANTI give consistent emissivities and it seems un- 
likely to us that both databases would give the wrong 
emissivities for such a relatively simple (He-like) ion. 
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Another possibility is that the underlying assump- 
tions for the calculation of the predicted line fluxes for 
FeXXV are incorrect. Our assumption of constant pres- 
sure logP = 13 implies that FeXXV is formed in an 
environment of logn e ~ 5 — 6, thus a rather low den- 
sity. However, in order to make a significant difference 
in the predicted FeXXV lines the density would have 
to be in excess of 10 16 cm -3 , and we reject this pos- 
sibility (see also lower part of Table [4} . The Fe XXV 
flux could be enhanced by resonant scattering into the 
line of sight if the plasma is not optically thin. This 
would only affect the resonance lines, and the FeXXVl 
line would also have to be enhanced but it is not detected. 
While the emissivities used to compute the Fexxv lo- 
cus curve uses the combined emissivities from the reso- 
nance, intercombination, and forbidden lines, contribu- 
tions from unresolved satellite lines are neglected. These 
can dominate the Fexxv complex at temperatures be- 
low the peak formation t empe rature (i.e., at logT < 
7.8; lOelgoetz fe Pradhanl l200lh . Lastly, a number of 
non-equilibrium pro cesses have significant effec ts on the 
Fexxv complex ( see lOelgoetz fe Prad han 2004). For ex- 
ample, recombination into excited states would enhance 
the FeXXV lines at the expense of Fexxvi and could 
explain the unusually high locus of the Fexxv lines. 
While significant effects of recombination of Fe xxvi into 
Fexxv violates the underlying assumption of collisional 
equilibrium, this does not necessarily imply that the 
other lines are also affected. Recombination affects only 
the ionization stages whose ionization energy is higher 
than the kinetic energy of the hottest plasma compo- 
nent. The energy required to ionize FeXXV into FeXXVl 
is 8.8 keV (equivalent to 10 s K), and that is higher than 
the hottest plasma component found from the APEC 
models of logT = 7.69 (Table |3J}. Meanwhile, the 
ionization temperature of Fexxiv into Fexxv is only 
2.04 keV (10 7 - 38 K) which is clearly lower than the hottest 
plasma component, and Fexxiv and Fexxv can be con- 
sidered to be in equilibrium. Also, S xvi and S xv are 
in equilibrium, since the ionization energy is 3.224 keV 
(= 10 757 K). We therefore conclude that only FeXXV 
and Fexxvi might not be in equilibrium while all other 
lines are, and our underlying assumptions are valid for 
these lines. 

In Table [HI we list the measured and predicted line 
fluxes for our best EMD for day 13.8 (Fig. [12]), giving 
element, rest wavelength, measured fluxes (after deblend- 
ing and correction for Ah — 5 x 10 21 cm -2 ), and ratios 
of measurements and predictions. We assume a value of 
Ah that is lower than that found from the APEC models 
(Table [5]) because we are unable to find an EMD model 
with that value that gives such good reproduction of all 
line fluxes. Also, with lower values of N-r we are hav- 
ing difficulties to find a good EMD model, and values 
as low as the interstellar value of Ah = 2.4 x 10 21 cm" 2 
can be excluded. We note that the EMD modeling is 
not an ideal way to determine Ah, and we refrain from 
determining a confidence range for Ah but consider it an 
uninteresting parameter, thus concentrating on the ele- 
mental abundances. We note that the chosen value of 
N-g is consi s tent w ith the total column density found by 
iBode et al.1 ([2006] ). 

In Table [7] we give the correction factors applied to the 
abundances. The EMD is scaled to reproduce the oxy- 



gen lines assuming solar O abundance (correction is 1.0), 
and the correction factors are thus equivalent to the re- 
spective abundances relative to solar. We also give the 
logarithmic abundances in the standard notation and list 
abundances relative to Fe (computed from the respective 
values relative to oxygen) in the bottom of Table [7] We 
determine the uncertainties of abundances by stepping 
through a grid of values, each time readjusting the EMD 
and computing a value of x 2 from the measured fluxes 
with their uncertainties for the selected lines. The listed 
uncertainties are derived from increases of \ 2 by one, 
which in the case of a 1-parameter model would be the 
1-<T uncertainties. We note that our model is not a pa- 
rameterized model. 

In the absence of any hydrogen lines, absolute abun- 
dances can only be determined via the strength of the 
continuum relative to the lines (see above). Since the 
mean EMD is not well constrained at high temperatures, 
the shape of a continuum model predicted by the EMD 
disagrees with the observed spectrum. It is not possi- 
ble to adjust the EMD without conflicts with some of 
the emission lines, and we thus refrain from determining 
absolute abundances with this method. 

We apply the same approach using the line fluxes from 
other observations listed in Table [3l and show the results 
in the bottom panel of Fig. [12] for day 111.7. All emis- 
sion measure values are significantly lower, and there is 
no indication for plasma hotter than logT = 7. The 
abundances are given in the last column of Table [6] and 
they are all consistent with those found on day 13.8. This 
conclusion also holds for the other datasets, although it 
is more difficult to distinguish between different EMD 
models. The reasons are lack of lines formed at low tem- 
peratures for the observations taken on days 26.1, 39.7, 
54.0, and 66.9 and lack of lines formed at high temper- 
atures for the observations taken after day 111.7. Also, 
some lines are blended with other nearby lines which can 
be disentangled with the Chandra HETGS, but not with 
the XMM-Newton RGS and Chandra LETGS. 

In Fig. [13] we show the minima of the emission mea- 
sure loci for all observed line fluxes using different plot 
symbols for each observation as explained in the legend. 
To guide the eye we connect the datapoints belonging 
to the same observations with solid and dotted lines in 
order that the evolution of the temperature structure 
can be identified. A sim ilar plot has been presented by 
iSchonrich fe Nesd (|2008h who assumed the same value 
of Ah = 2.4 x 10 cm' 2 for all observations and solar 
abundances. They found significantly higher loci for the 
N and O lines for days 39.7 and 66.9 compared to all 
other observations. Since during the SSS phase these 
lines appeared on top of the SSS continuum, they con- 
cluded that these lines are formed within the outflow or 
are at least affected by the SSS radiation. We now use the 
new abundances from Table [7] and Ah = 5 x 10 21 cm" 2 
for days 13.8 and 26.1 and N K = 2.4 x 10 21 cm" 2 for the 
rest (see Tabled]). The higher value of Ah for days 13.8 
and 26.1 leads to higher loci of the low-temperature lines 
of O and N which are now consistent with those mea- 
sured for day 39.7 and 66.9. In order to attribute these 
lines to either the outflow or the shock therefore requires 
accurate knowledge of the value of Ah. But even with 
our improved measurements of Ah the situation remains 
ambiguous because, if the O and N lines observed on top 
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of the SSS continuum are formed somewhere within the 
outflow, they might be subject to higher values of iVn as 
is suggestive from the shock plasma. This would increase 
the discrepancies again, and with these uncertainties, it 
may never be possible to decide whether these emission 
lines are formed in the outflow or in the shock. 

At the other end of the temperature distribution, at 
log T > 6.8, a steady decrease of emission measure can be 
identified. The slope of the temperature distribution to- 
ward the highest temperatures seems to become steeper, 
indicating that the hot plasma is cooling rapidly, while 
below logT = 6.8 the line emission measures also de- 
crease, but the slope remains about the same, and the 
cool component thus cools at a slower rate. 

4.4. Comparison of models 

While the APEC models introduced in §4~2l include all 
available atomic information, the line-based approach 
used in £14.31 has the advantage that only the most 
reliable information is selected. Less certain lines 
(e.g., lines with transitions involving higher principal 
quantum numbers) are discarded. In £14.31 the crude 
assumption of constant pressure is sufficient because 
only lines that are not density-sensitive are selected. 
The APEC model makes the even cruder assumption o f 
a low-density plasma (logn e = 0: iSmith et al.l l2001ah . 
while density-sensitive lines are not excluded. However, 
most density-sensitive lines are weak. Another difference 
is the assessment of the goodness of a model. In our 
line-based approach the fitting yields best reproduction 
of measured line fluxes (taking line broadening into 
account), while with XSPEC the count rate in each spec- 
tral bin has to be reproduced. In cases where lines are 
present in the observed spectra that are missing in the 
atomic database, we can ignore them with the line-based 
approach, while with XSPEC the existing lines in the 
atomic database can be used to force an acceptable fit 
of these spectral regions. Furthermore, the individual 
emission lines contribute only to a few spectral bins 
such that negligence of emission lines is badly penalized 
compared to negligence of the continuum. On the other 
hand, the APEC models are much better suited to assess 
the hottest temperature component via the continuum, 
and the corresponding shock velocities can only be 
derived from the APEC models (see Fig. [TTj). For more 
discussion on line-base d and global fitting approaches 
we refer to iNessI (|2006T ) . We thus need both approaches 
for robust conclusions. 

In the top panel of Fig. [12] we include the temperatures 
and volume emission measures of the three components 
derived from the 3-T APEC models fitted to the spec- 
tra of day 13.8 for comparison with the mean EMD. All 
three components yield higher values of emission measure 
than the continuous EMD, because all emission from the 
smooth distribution of the EMD model is concentrated 
in only three isothermal components. The emission mea- 
sure of the hottest component of the APEC model is 
much higher than in the EMD model. Since not enough 
emission lines are formed above 10 7 ' 5 K to constrain the 
EMD, the line-based approach is clearly inferior in this 
temperature regime, and the hottest component is driven 
by the continuum. We further find that only the lines 
that are formed at high temperatures are affected by re- 



combination (see SI4.3P . The second temperature compo- 
nent of the APEC model coincides with the temperature 
of the hotter peak of the mean EMD. The third compo- 
nent is closer to the minimum between the two peaks of 
the mean EMD (see Fig. [T2|) . The higher temperature in 
the APEC model could explain why the N abundance is 
higher than that derived from the EMD modeling. Since 
each temperature component is isothermal, the N lines 
are not as efficiently produced by the cool component, 
which has to be compensated by a higher N abundance 
in order to fit the N lines. We also note that the Nvi 
line at 28.78 A, which is clearly detected (see Tabled]), is 
completely ignored by the APEC models but represents 
an important constraint on the EMD model. Finally, in 
the XSPEC fits we could not account for line broadening 
(see §4.2]) . while we used the line fluxes integrated over 
the entire profile for the EMD reconstruction method. 
We therefore regard the abundances of O and N derived 
from the EMD model as more reliable than the values 
derived from the APEC model. 

For all models, we give only statistical uncertainties. 
Uncertainties from the atomic physics are not included 
and are a source of additional systematic uncertainty 
( H4.5p . Since many lines with poorly-known atomic 
physics are included in the APEC models the system- 
atic uncertainties of the APEC models are higher than 
those from the EMD models. 

For comparison of the elemental abundances derived 
from the two approaches, we rescale those obtained 
in §4.21 because the referen ce abundances used in th e 
APEC models are those by lAnders fe Grevessd (p.989), 
while for the EMD reconstructi o n me thod we have 
used those by iGrevesse fe Sauval (jl998). We rescale 
by a factor (Ax/A r o) g rev/(A r x/A r o)and, where the sub- 
script s 'grev' and 'and ' deno te th e abundance ratios 
from IGrevesse fe Sauvall <| 19981 ) and lAnders fe Grevessel 
(1989), respectively. The correction factors are 0.933, 
1.23, 1.259, 1.259, 1.66, and 0.85 for N, Ne, Mg, Si, S, 
and Fe, respectively. In Table [7] all derived abundances 
are listed for comparison. While the abundance ratios 
relative to O are discrepant, those relative to Fe agree 
much better, except for N/Fe. We attribute these dif- 
ferences to the low formation temperatures of the N and 
O lines. We attribute these differences to the less cer- 
tain N and O abundances derived from the APEC mod- 
els that underestimate the amount of cool plasma (see 
above). We expect overestimated N and O abundances 
and possibly also Ne in the APEC model which explains 
that all abundances relative to O are lower in the APEC 
model compared to the EMD model. Since the N lines 
are formed at lower temperatures than the O lines, the N 
abundance is affected to a higher degree. The Fe lines are 
formed over a large range of temperatures and are there- 
fore not as strongly affected, leading to the overall better 
agreement of all ratios relative to Fe. For the observation 
taken on day 111.7 only the EMD modeling yields some 
constraints on the elemental abundances, which demon- 
strates the strength of this approach over the spectral 
fitting. 

4.5. Uncertainties 

While the results from §2] are directly based on the 
observations, all results from this section, §3 depend on 
model assumptions which are described in £14.11 All error 
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Fig. 13. — Volume emission measures (VEM) at the peak formation temperature for selected strong lines, assuming the abundances 
listed in Table [7] The connecting lines guide the eye, but are not EMD models. The plot symbols indicate the time of observation after 
outburst according to the key in the bottom right legend. We use dark gray for observations taken before the SSS appeared, light gray for 
observations during the SSS phase and black for after the SSS spectrum had disappeared. 



estimates given in this paper are statistical 1-ct uncertain- 
ties which give the 68.3-per cent probability that fitting 
the same model to a new observation with the same in- 
strumental setup results in parameters within the given 
uncertainty ranges. They thus only describe the precision 
of our measurements, but not the accuracy (sum of sta- 
tistical and systematic uncertainties) which depends on 
the calibration of the observations but also on the choice 
of a model. For comparisons of measurements taken with 
the same instrument (given in Table [l}, the systematic 
calibration uncertainties can be neglected but have to be 
kept in mind for absolute numbers and comparisons with 
different instruments (e.g. the flux evolution shown in 
Fig.[|5]). For line profiles (Fig. [7]) and line ratios (Fig. [5]), 
the cross-instrument calibration uncertainties are negli- 
gible. 

For the flux measurements presented in Fig.[SJ system- 
atic uncertainties arise from the choice of band width 
which excludes a fraction of the total X-ray emission. 
Owing to the presence of the SSS spectrum between days 
39.7 and 66.9, the contribution from emission between 11 
and 38 A to the shock emission can not be determined, 
yet the fraction of soft emission may be higher compared 
to before day 39.7. Attempts to determine this contribu- 
tion from the XSPEC models failed because we have no 
constraints from observations because the much stronger 
SSS emission dominates at long wavelengths. We esti- 
mate that the systematic uncertainties on the Swift light 
curve shown in Fig. [5] are small, because the spectral 
shape hardly changes after day 100 (see Table[5]), and the 
X-ray flux thus scales directly with the observed count 
rate. We note that systematic uncertainties from direct 
rescaling are smaller than the method of flux determina- 
tion via model fitting to each individual Swift spectrum. 

For the line shift measurements presented in Fig. [7] 
and i j3.ll systematic uncertainties can arise from line 
blends and background noise, which affects the weaker 
lines more than stronger lines. For the measurement of 
line fluxes we have not accounted for uncertainties from 
choice of a source background (c ontinuum) underneath 
the lines and from the line widths. iNess fe Jordan! (|2008f ) 



found that the uncertainties in the line widths have less 
effect than the choice of continuum. For the line ratios 
presented in Fig. [8l additional systematic uncertainties 
from the continuum are less than 5 per cent if the contin- 
uum is assumed to be uncertain at the 20 per cent level. 

For the model fits, the main sources of systematic un- 
certainty are the model assumptions and uncertainties 
in the atomic data. Both are difficult to quantify. The 
hottest temperature components of the 3-T fits presented 
in Table [5] are driven by the continuum. The Fe xxv 
and FeXXVl lines are the only strong lines formed at 
temperatures higher than ~ 10 7 K, and their absence 
in a model poses little resistance against an overesti- 
mated continuum temperature. For observations with 
weak continuum, arbitrarily high temperatures can thus 
result, and we give higher confidence to the temperatures 
of the hottest component determined from observations 
before day 39.7 (see Table [5]). In view of these large sys- 
tematic uncertainties of the later observations, the trends 
of shock velocities plotted in Fig. [TT] are only reliable for 
the early observations. Since the apparent f -1 / 3 trend 
that is suggestive from the values derived for day 39.7 
and 111.7 is based on the higher temperature estimates 
compared to the expected i -1 / 2 trend, the conclusion of 
adiabatic cooling instead of the theoretically expected 
radiative cooling has to be treated with caution. 

The uncertainties in elemental abundances result- 
ing from the emission measure modeling (Table [7]) 
are extremely difficult to quantify. They mainly de- 
pend on the uncertainties in the line flux measure- 
ments and the uncertainties in the determination of 
the mean EMD. INess fc Jordan] (2008) found from com- 
pa rison to independent analy ses of the same spectrum 
by iSanz-Forcada et"aLl (|2003f ) that the abundances are 
quite robust against changes in the assumed mean EMD. 
Since the abundance determination is based on strong 
lines originating from few-electron ions (thus well-known 
atomic physics), we are confident that the given uncer- 
tainties that are based on the measurement uncertainties 
of line fluxes are realistic. 
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5. DISCUSSION 

The X-ray grating spectra give the deepest insight into 
the properties of the shocked plasma. With 12 grating 
observations we can follow the changes in the X-ray flux 
and plasma temperature, and we have determined the el- 
emental abundances from the emission lines. The X-ray 
flux can directly be integrated over the spectrum without 
the need of a model and can thus be considered an ob- 
served quantity. We are further able to determine tem- 
peratures independently from line flux ratios and from 
spectral models, and we are able to detect different tem- 
perature components. The ability to measure line fluxes 
also enables us to pursue two independent approaches to 
determine the elemental abundances, yielding robust re- 
sults. We confirm that the underlying assumptions are 
generally satisfied, but we find evidence for recombina- 
tion in the Fexxv He-like triplet lines and the Fexxvi 
lines. Since these are the only ions whose ionization en- 
ergy exceeds the kinetic energy of the hottest plasma 
component, we regard the assumptions discussed in $4] 
valid for all other lines. 

We observe a power-law decay in X-ray flux, and the 
power-law index changes from a = —5/3 before day 
~ 70 to a = —8/3 during the later evolution (see 
Fig. [5]). While the t~ 5 ^ 3 decay is predicted for a ra- 
diatively cooling shock traversing an r~ 2 density distri- 
bution, the change to a t -8 / 3 decay at late times cannot 
be sim ply explained. However, note that IVavtet et al.l 
(|2007h found that in some circumstances the forward 
shock could decelerate faster than the standard cooled 
momentum-conserving models predict. In addition, it is 
interesting to note that a steeper decay may result from 
the breakout of the forward shock into a lower density en- 
vironment when it reaches the edge of the red giant wind 
re-established in the 21 year interval between explosions. 
Indeed, simple calculations suggest that this breakout 
may happen around t he time of the apparent break in 
decay curves (see, e.g. jMason et a,l.lll987t lO'Brien et all 
1992). One possible problem with this conclusion how- 
ever is that the immediate post-shock temperature in- 
creases, rapid adiabatic expansion of the emitting mate- 
rial will likely lead to a drop in the temperature derived 
from fits to the X-ray spec tra, which does no t seem to 
be the case. Furthermore, iBode et all (|2008f ) conclude 
from the apparent early onset of the adiabatic expan- 
sion phase in the forward shock that shock breakout will 
occur much later than realized from X-ray observations 
that in 1985 started only at 55 days post-outburst. 

The two peaks in the emission measure distribution 
(see Fig. [12]) could represent the forward and reverse 
shocks, where the cooler component wo uld be the re- 
verse shock. But according to models bv lBode fc KahrJ 
(198$, the reverse shock may not be a significant source 
of X-ray emission at this time. Alternatively, the two 
temperature regimes could relate to forward shocks prop- 
agating into two distinct circumstellar density regimes, 
e.g., denser equatorial regions and less dense polar re- 
gions as sugg ested by VLBA/EVN and HST observations 
(jBode et al.M2007t lO'Brien et all 120061) . More realistic 
hyd rodynamic modellin g is needed, and this is underway 
(see IVavtet et ai1l2007D . 

In several Classical Novae, e.g., V382Vel, X-ray spec- 
tra similar to those of the late evolution in RS Oph 



have been obtained after the nova had turned off and 
were attributed to t he ejecta that are radiatively cool- 
ing (|Ness et al.ll2005D . It is thus possible that the cooler 
component observed after day 70 originates from the ex- 
panding ejecta, while the early hot and cool components 
represent the shock emission. 

5.1. Elemental Abundances 

For our EMD models and multi-temperature fits, we 
have assumed uniform abundances throughout the X-ray 
emitting plasma. If the relative contributions from each 
component to the observed X-ray emission change with 
time and their compositions are significantly different, 
then in principle, changes in observed abundances de- 
rived from different observations are possible. However, 
the relative abundances obtained for days 13.8 and 111.7 
do not differ significantly, and we conclude that no signs 
of changes in the composition are detectable. We there- 
fore focus on our measurements for the best dataset taken 
on day 13.8. There are several issues for which the accu- 
rate determination of elemental abundances is of interest. 
The composition of the WD can be used as an argument 
whether or not RS Oph will explode as a canonical SN la. 
The WD may be near the Chandrasekhar limit, because 
of the short recurrence time scale, but if the WD con- 
sists of too many heavy nuclei (Mg, Ne rather than C, 
O), the amount of available nuclear binding energy is not 
sufficient to totally disrupt the WD. Note that it is the 
explos ions of CO WDs t hat best fit the observations o f 
SNIa (|Leibundgutil2000tlHillebrandt fc Niemeverll2000ft . 

Since hydrogen is fully ionized in X-ray emitting 
plasma no H lines can be observed, and we cannot de- 
termine absolute abundances. In both of our approaches 
the attempt to determine absolute abundances from the 
strength of the continuum relative to the lines leads to 
ambiguous results, and we thus focus on relative abun- 
dances. 
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Fig. 14. — Photospheric abunda nce ratios of 20 M giants given 
by ISmith fc Lambertl (1985, 1986). The star symbol is our mea- 
surement for RS Oph as listed in Table [7] Uncertainties are given 
by gray error bars. 



GNO -cycled material 

There is no doubt, from both of our approaches, 
that nitrogen is overabundant. This points to CNO- 
cycled material, which could have been produced dur- 
ing the outburst. However, the companion is a red 
giant and is as such li kely overabundant in ni t rogen 
as well (see Fig [H a nd ISmith fc Lambertl 1 19851 119861 : 
iSchmitt fc Ness! I2002D . The acc reted material on the 
WD is thus already N-enhanced (|Sniiderslll987af ). as is 
the wind into which the ejecta run. In addition, nu- 
clear burning during the outburst can contribute with 
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more CNO-cycled material. Furthermore, N-enhanced 
material from previous outbursts may al so contribute 
to the X-ray emitting plasma. Note that IWalder et al.1 
(2008) have determined from theoretical models that the 
ejecta are dominated by the RG material. In Fig. [14] 
we show the N a bundance measurements of 20 M gi- 
ants presented by ISmith fc Lambertl (|1985l . 119861) ver- 
sus their O measurements. The star symbol indicates 
our result for RS Oph, and it is clearly higher than 
any of the M giants, while the O abundance is consis- 
tent with most of the M stars. The comparison with 
other M giants thus suggests that CNO-processed mate- 
rial may have been added by the outburst. In support 
of this, the N/O abundance ratio derived from our two 
ap proaches is roughl y simi lar to the value determined 
bv ISchonrich fc Nessl p008T ) ([N/O] = 0.7), who deter- 
mined the ratio from the line column densities in the 
SSS spectra. Our result is also consistent wi th the abun- 
dance measurements bv lContini et al.l (|1995h which yield 
an abundance ratio of [N/Fc] = 1.18. The photospheric 
abundan ce measurement s of th e RS Oph secondary pre- 
sented bv lPavlenko et all ([20081 ) yield [N/Fe] = 0.9 which 
can be compared to our abundance ratio [N/Fe] = 1.14. 
Th is is lower than ei t her ou r value or that determined 
by iSchdnrich fclfcsi (120081 ) ([N/Fe] = 1.0 under the 
assumption of [O/Fe] = 0.3), suggesting that the out- 
burst has added about 20-40% of nitrogen into the X-ray 
detected mater i al. W e note that [N/Fe] as derived by 
iPavlenko et al.l (|2008f ) is higher than any of the values 
plotted in Fig. [TU This indicates that the secondary in 
the RS Oph system is not a typical M giant which may 
be caused by the fact that the RS Oph secondary has lost 
significant amounts of material from its outer envelope 
during accretion to the WD. 

If the material was CNO cycled, then we would 
expect it to be und erabundant in carbon (see, e.g., 
ISchmitt fc Nessl 120021 ). but all carbon lines are at wave- 
lengths where strong effects by interstellar and circum- 
stellar absorption dominate. Moreover, uncertainties in 
Nn propagate to large uncertainties in any line fluxes or 
upper limits of the C lines (see Table [3]). For any quan- 
titative assessmen t of CNO b urning, the carbon abun- 
d ance is essent i al ([D rake 20031). 

iNelson et al.l (|2008f) discussed the carbon abundances 
and obtained an underabundance of C/N by a factor of 
0.001 in order to explain the absence of the Cvi K-shell 
absorption edge in the SSS spectra at 25.37 A. This 
is much lower than the C/N abundance ratio of 0.05 
that one deri ves from the [C ] = —0 .4 and [N] = +0.9, 
estimated by IPavlenko et alj ([2008) for the secondary. 
Strong ca rbon lines were a lso observed in the IUE 
spectra bv IShore et al.l (|1996f) . The absence of the Cvi 
edge is not sufficient evidence for a C underabundance 
because other high-ionization absorption edges are also 
absent in the SSS spectra of RS Oph. For example, the 
ionization edges of Ovin at 14.23 A, Ovu at 16.77 A, 
Nvn at 18.50 A, and Nvi at 22.457 A are not present 
although these elements are not underabundant (see 
figs. 3 and 4 in iNess et al.l 120071 and figs. 8 and 9 in 
INelson et all 120081 ) . Since iNess et ail (|2007l ) measured 
considerable blue shifts in the absorption lines observed 
during the SSS phase of RS Oph, the shell around the 
WD must b e expa nding with high velocities. While 
INelson et al.l ([2008) found no blue-shifts in the same 



lines, they quote high velocities derived from line shifts 
of emission lines of C VI Ly series lines, however, without 
the presence of the Lya line (Table [3]). For an envi- 
ronment with high expansion velocities, the absorption 
edges may be washed out by the expansion. We therefore 
argue that without measurements of emission line fluxes 
of carbon, the C abundance cannot be determined from 
the X-ray spectra. Finally, a C/N abundance of 0.001 is 
far too low for any material that has undergone CNO 
cycle processing. The large N abundance r equires some 
processing back to C dStarrfield et al.ll200l . 

a elements (Z < 22) 

The abundances of the a elements Ne, Mg, Si, S are sig- 
nificantly higher than O and Fe using both approaches. 
We use O as the reference element for our analysis, and 
our Ne/O abundance ra tio is consist e nt the value found 
from IR observations bv I Evans et all (|2007bD . This indi- 
cates that either Ne is ove rabundant or O is underabun- 
dant. iContini et al.l (|1995f ) found significant underabun- 
dance of O/H and of Ne/H of - 10% solar. We also 
list the ratios relat ive to Fe in T a ble [71 The abundance 
determinations by IContini et al.1 (|1995f ) also yield high 
abundance ratios of [Mg/Fe] = 0.73 and [Si/Fe] = 0.86, 
which are roughly consistent with our results. Since Fe is 
neither produced nor destroyed in the outburst, we argue 
that all a elements are overabundant in the X-ray emit- 
ting plasma. Since we determined these values with two 
independent approaches using spectra with strong lines 
of relatively simple ions with well-known atomic data, we 
regard these results as reliable. 

There are three possible explanations: (1) the a ele- 
ments were produced during the outburst, (2), the com- 
position of the secondary, which provides the accreted 
material and dominates the composition of the stellar 
wind, is overabundant in these elements (3) WD core 
material is enriched in the a elements, and that material 
has been mixed into the outflowing material. 

(1) The production of a elements via nuclear burning 
requires high temperatures and densities. While the re- 
quired temperatures are likely reached during the TNR, 
the densities are too low in order to produce significant 
amounts of a elements, even if the WD mass is near the 
Chandrasekhar mass. 
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Fig. 15. — Photosphe ric abundance ratios (r elative to solar b y 
IGrevesse fc Sauva 1 1998) of 17 M giants given bv lRich et all (2007). 
The star symbol is our measurement for RS Oph as listed in Table[7] 
Uncertainties are given by gray error bars. 

(2) Direct abundance measurements of a elements 
for th e secondary are not available, and IPavlenko et ail 
(2008) only determined Fe, C, and N abundances. 
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iRich et all (|2007f) measured photospheric abundances for 
17 M giants in the inner bulge of the galaxy from IR spec- 
troscopy. In Fig. [15] we show their abundances of [Si/Fe] 
versus [O/Fe] (top) and versus [Mg/Fe] (bottom), rela- 
tive to lGrevesse fc Sauvall (jl998). All M giants are over- 
abundant in Si, O and Mg. Stars with a higher Si abun- 
dance are also higher in Mg and O, indicating different 
concentrations of the ashes of He shell flashes, proba- 
bly due to different stellar ages. Our Si and Mg abun- 
dances for RS Oph are significantly higher than any of 
the M giants in the sample, but the O abundance is not 
higher. The bottom panel of Fig. [15] shows that our value 
is significantly higher in both elements than any of the 
M giants in the sample, and the concentration of all a- 
elements could be higher. Since the outer envelope of the 
RS Oph secondary has been stripped away and accreted 
by the WD, higher abundances of all elements produced 
by He-burning can be expected. The surface composi- 
tion may resemble that of the inner regions of normal 
M giants. Meanwhile, the O abundance could reflect a 
balance between a higher O abundance produced by He 
burning and a lower O abundance in material that has 
a higher concentration of CNO processed material as ev- 
ident from Fig. Q3] Our measurements can thus reflect 
the composition of the accreted material. 

(3) If direct abundance measurements of the a ele- 
ments of the secondary are different from our values, 
then the only remaining possibility would be that the 
WD contributes significantly to the composition of the 
observed plasma. The composition of the WD is domi- 
nated by elements that have been produced in the pro- 
genitor star. In order to be observable, WD material has 
to be dredged up and mixed into the outflowing mate- 
rial. While this is possible, no observations have so far 
revealed any signatures of WD material in the emitting 
regions of RS Oph. 

5.2. Is RS Oph a SN la progenitor? 

There is some speculation that RS Oph is a SN la pro- 
genitor. The short recurrence time scale suggests that 
the underlying WD is of h igh mass, possibly near the 
Chandrasekhar limit (e.g.. iDobrzvcka fc Kenvonl 11994 
IShore et all 119961 : iFekel et alj I2000D . lHachisu k Katol 
42001D reported a high-mass WD of ~ 1.35 ± 0.01 M 
(for both, high- and low-metallicity models), but their 
modelling neglects the environment and is not real- 
istic. An eve n high er W D mass was e s timate d by 
ISokoloski et~aT1 (|2006l ). and lHachisu et all (|2007l ) also 
found that the WD is growing in mass. According to 
their estimated growth rate, the Chandrasekhar limit 
would be reached in a few times 10 5 years. 

ONe WDs do not provide enough nuclear binding 
energy (see, e . g., to tal binding energy ca lculations by 
iGamezo et all 120031 or ICalder et all I2007I) . and when 
reaching the Chandrasekhar mass limit, they implode in 
a core collapse without an explosion. The explosion mod- 
els of SNIa predict that the WD has to be a CO WD for 
the supernova to be a canonical la explosion. However, 
the WD composition is not unambiguously determined. 

A strong argument against RS Oph being a SN la pro- 
genitor is th e large amount of hy drogen and helium in the 
system (e.g. lEvans et al.ll2007aJ ) that has to be removed 
before the WD reaches the Chandrasekhar limit in order 
for the SN la outburst to be consistent with observations 



of hydrogen-deficiency (Filippenkol ll997l ) . While the sec- 
ondary loses hydrogen during the evolution of accretion 
and repeated outbursts, it would have to be considered 
too much of a coincidence if all the hydrogen is consumed 
in nova outbursts at the same time as the Chandrasekhar 
limit is reached. Considering these points, we believe 
that it is unlikely that RS Oph is a SNIa progenitor. 

6. SUMMARY AND CONCLUSIONS 

The high-resolution X-ray spectra of the 2006 outburst 
of RS Oph provide unique insights into the properties and 
evolution of the outburst. The properties are best deter- 
mined from the dataset taken on day 13.8. At this time 
of the evolution the nova was brightest in hard X-rays, 
and we have the best coverage with almost simultaneous 
Chandra and XMM-Newton observations. From these 
observations we have derived the elemental abundances 
with an overabundance of N, Ne, Mg, Si, and S, relative 
to Fe, indicating that this material has undergone CNO 
burning and He burning, respectively. Both processes 
can occur in the RG companion star, and our observa- 
tions could reflect the donor material. CNO burning also 
occurs during the outburst, and by comparison of our re- 
sults with direct measurements of the N abundance for 
the secondary bv lPavlenko et all (|2008l ). we estimate that 
about 20-40% of the nitrogen could have been produced 
during the outburst. Meanwhile, no a elements can be 
produced during the outburst, but the underlying WD 
might be enriched in a elements. 

We find Mg and Si si gnificantly h igher than in any M 
giant in the sample bv lRich et alJ ((2007) which can ei- 
ther be explained by the RG being different from other M 
giants, or WD material has been dredged up during the 
outburst and mixed int o the e.jecta. Similar ly high val- 
ues have been found by IContini et alj (|1995f ) . Since the 
secondary has likely been stripped off it's outer layers 
by mass loss onto the WD which has then been ejected 
during the outburst, the composition of the ejecta may 
be sampling deeper layers of the M giant companion as 
compared to the photospheric composition of normal M 
giants. In order to determine whether WD material is ob- 
served, direct measurements of the Mg and Si abundance 
in the photosphere of RG are needed for comparison. 

We also determined the evolution of the temperature 
structure. The early observations show both a hot and a 
relatively cooler plasma component, while the later ob- 
servations only display a cool component. Both com- 
ponents decay with time. Before day ~ 70, while the 
plasma is dominated by the hot component, the decay 
rate is slower than later in the evolution. The early de- 
cay rate is consistent with a radiatively cooling plasma 
while the later evolution can be explained by the expan- 
sion of the ejecta. The temperature evolution derived 
from models are consistent with radiative cooling of the 
hot component. 

While it has been suggested that the WD is close to 
the Chandrasekhar limit, we note that there is currently 
too much hydrogen in the system for a supernova explo- 
sion to satisfy the spectroscopic features of a la explosion 
(Filippenko 1997). The hydrogen and helium will have to 
be removed before the Chandrasekhar limit is reached. 
Further, the underlying WD has to be a CO WD, but 
the composition of the WD is difficult to determine since 
WD material has to be dredged up and ejected. Our 



23 



abundance analysis suggests that this may be the case, 
but the composition of the accreted material is not well- 
enough known for solid conclusions at this time. 
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APPENDIX 
Statistics 

Some of the spectra are sufficientl y bright for s tandard % 2 minimization, while others are extreme ly faint, and low- 
count statistics, as recommended by ICashl (|1979() . have to be applied. According to Wilks theorem (|Wilkslll938l ). the 
Maximum Likelihood (ML) technique converges to the \ 2 statistics in the high-count limit, and we thus apply the ML 
technique to all observations. The likelihood is defined as 

N 

C = — 2^~"^(rtj km, — m,) (f) 

i=i 

with rii being an observation and m, a model, both defined on the same grid of N bins. According to Wilks theorem C 
serves as a goodness criterion in an equivalent way as x 2 , but it is derived from the Poissonian probability distribution 
function (PDF) rather than a Gaussian PDF. In order to conserve the Poissonian nature of the data, the instrumental 
background must not be subtracted. Instead, the model has to be added to the background for comparison with the 
raw data, because the sum of two Poissonian statistics is Poissonian, while the difference is not (jCashlll979| ). 



To assess the statistical uncertainty ranges of each parameter in a multi-parameter model we compute the curva- 
tures of the likelihood at the respective best-fit values from the second derivatives. In order to assess the correlated 
uncertainties, we also compute mixed derivatives, thus the full Hessian matrix which, for n free parameters, is an n x n 
matrix with the partial second derivatives of the likelihood curve, C = C(Ai, Aj), 



dC 2 
dAidAj 



Hi i = (2) 



where Ai and Aj represent two free parameters with i and j ranging from 1 to n. We determine Eigenvectors, EVi(Aj) 
and Eigenvalues, ev(Aj), and compute the correlated uncertainties from 



3 mr/ 1 12 



J ev(Aj) 

i— 1 v J ' 

with AC = 2.3 and AC = 3.53 for 1-er uncertainties in the cases of two and three free parameters, respectively ()Stronel 
1985). Uncorrclatcd errors can be computed by setting all off-diagonal elements to zero before computing Eigenvectors 
and -values. 



